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ABSIHAOT 


In this investigation the transfer matrix approach 
is used to determine the natural freq.uencies and mode 
shapes of a rotating large aspect ratio asymmetric 
tapered turbine blade. Some of the predominant factors 
like shear deformation, rotary inertia, root flexibility 
and stagger angle which affect the natural freq.uencies 
of the blades are also studied. 

Hunge-Eutta method is used to solve the matrix 
of the variable co-eff ecient s of system of- first-order 
differential equations of motion for coupled bending- 
torsion and pure flexural mod.es.' 

Typical results are presented in the tabular as 


well as graphical form 



CHAPTER I 


IKTRODUCTIOH 

1.1 DescriiPtion of the Problem 

Knowledge of the resonant characteristics of 
turbine and compressor blades is of utmost importance 
in insuring reliable long-life- turbine engines. The 
damage caused by blades failing due to vibratory 
fatigue can be catastrophic at worst and at the very 
least result in additional engine development costs due 
to redesign and repair. Comprehensive analytical and 
experimental methods for analyzing the vibrational 
characteristics of turbine blading are essential in the 
avoidance and prediction of turbine engine blade vibration 
problems , 

With the large increase in the use of axial flow 
compressor 'and turbine, determination of natural frequencies 
and modes of vibration is of considerable importance in 
the design of blades used in turbomachines. It is known 
from earlier investigations that most of the blade failures 
occur due to fatigue at the blade roots, and it occurs 
when vibrations take ^ at or near resonant conditions. Eor 
the safe design, the designer has to make sure that the 
natural frequencies obtained by theoretical calculations 
are reasonably away from the forcing frequencies. 

The difficulties facing the design of blad^g are 
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many. The magnitude of the problem lies in the fact that 
in some turbomachines there can be as many as a few 
thousand fixed and rotating blades of different characteristics 
and the failure of even one of them will force a shutdown. 

The major problem of design is a fair prediction of the 
natural frequencies of the blading at the earlier stages 
by theoretical calculations. 

In* determining the natural frequencies, the designer 
can idealize the large aspect ratio blade as a uniform 
cantilever beam and find the effects of the following 
primary important factors : 

1, Rotary inertia and Shear deformation, 2. Rotation of 
the disk, 3* Taper of the blade cross-section, 4. Staler 
angle of the blade, 5« Root flexibility, 6, Disk elasticity, 

7. Disk radius, 8, Pretwist of the blade, 9* Asymmetry of 
the blade, 10. High temperature, 11. Aerodynamic excitation, 

12. Aspect ratio of blade. 

As turbomachine blades are often joined together , either 
at their tips — shrouding or at some intermediate location 
-- lacing. In addition to the above mentioned parameters 
few more parameters that affect the response of blades 
are as follows : 

1. Shrouding, 2. Lacing wires, 3. Large amplitude, 4. Damping, 
5. Material. 

A good deal of theoretical and experimental work is 
reported in the literature which enables the designer to • 
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calculate the natural frequencies, but it is still a long 
way to go before the last word can be said on this problem* 

As one cannot entirely depend upon the theoretical 
calculations, it becomes essential to mate full~scale 
model tests on the prototype and study the vibration 
characteristics, i.e., the natural frequencies euid their 
corresponding nodal patterns. But the experimental 
procedure is time-consuming and it is very expensive to 
conduct tests on a number of designs. Thus it becomes 
essential for the disigner to have theoreMcal knowledge 
as well as practical means for conducting model tests in 
order to .predict the vibration characteristics of the 
blading at an earlier stage of design. 

The most important modes of vibration from practical 
considerations, are known to be the first three flapwise 
bending modes, the first chordwise bending mode and the 
first two torsional modes. Depending on the characteristics, 
the turbine blading can execute either uncoupled bending 
or torsional vibrations or coupled bending torsion vibrations. 
Coupled bending torsion vibrations occur when the centre 
of flexure does not coincide with the centroid of the 
blade cross-section. When the blade is pre twisted, vibrations 
are further coupled between the two bending modes. The 
differential equations of motion for the case of a pretwisted 
blade with assrmmetric cross-section are quite complex and 
closed form solution of such equations is difficult to obtain. 



Also the natural frequencies determined are affected by- 
shear deformations, rotary inertia, warping of the cross-£. 
section, root fixing and coriolis accelerations. In 
general the equations of motion will be six coupled partial 
differential equations coupled between the two bending 
deflections, the two shearing deflections, the torsional, 
deflection and the longitudinal deflection. Further the 
warping function of the blade cross-section will be obtained 
by a modified Poission equation, taking into account the 
dynamic conditions of the blade. Since tfee longitudinal 
natural frequency of the blade is generally very high 
and there is less possibility of exciting a lofigitudinal 
mode of the blade and hence this deflection may be ignored. 
Carnegie [1] considered the coriolis accelerations due to 
the inward displacement of the blade element and derived 
the nonlinear partial differential equations which are very 
complex, it would be a futile exercise in the present day 
to attempt a solution of them taking into account all the 
effects. When the blades are shrouded, the coupling, 
between the blade and the shroud is to be further taken 
into account. Thus, theoretically it is an uphill task 
to determine the natural frequencies of an actual turbine 
blade with all the effects mentioned. 

As stated above, the number of parameters which affect 
the response of the blades are quite large and is rather 
impossible to consider all these aspects at one time. 
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Fxirthermore, to evaluate their relative importance, it 
is easy to incorporate these parametersone hy one* 

In this investigation, transfer matrix approach is 
used to determine the natural frequencies and mode shapes 
of a rotating hladf on the low temperature side of the 
turbine — idealized as a large aspect ratio cantilever 
beam tapered in both depth and width with different support 
flexibilities. The formulation of the problem is derived 
for pure bending and coupled bending torsion modes of 
vibration — including some of the predominant factors 
like shear deformation, rotary inertia, root flexibility, 
stagger angle etc* 

Chapter II contains up- to-date literature survey 
on most of the parameters considered to be effective for 
affecting the response. Chapter III and 17 cover the 
formulation of the problem for pure bending and coupled bending 
torsion modes of vibration respectively. Numerical results 
and conclusions are presented in Chapter 7. Chapter 71 
deals with the recommendations for further study in this 


direction. 



CHAPTER II 


2.0 LITERATURE SURVEY 

In the last twenty five years a large number of 
technical papers have appeared, describing the different 
approaches for the vibration analysis of turbomachinery 
blades. The general methods of solving the blade vibration 
problems are : 

1. Exact method, 2. Raylei^' s method, 3. Rayleigh- Ritz 
method, 4. Ritz and Gralerkin method, 5. Matrix methods, 

6. Prohl method, 7. Numerical methods to determine the 
frequency equation, 8, Integral equation approach, 9* Finite 
difference method, 10. Finite element method, 11. Perturbation 
method, 12. Frequency response methods* 

The classical approach of solving the differential 
equations of motion of a cantilever blade is discussed in 
most of the books on vibration, -The solution is possible only 
under simplifying conditions. Matrix methods are considered 
by Duncan [2], Duncan and Collar [3], Dokumaci [4], Carnegie, 
Thomas and Dokumaci [5] and Scanlan and Rosenbaum [6], 

Prohl [ 7 ] method is very successful for bending vibration 
problems of the blades. 

Based on Holzer and Myklestad equations, Rao [8,9] 
developed numerical procedures to determine frequency equations 
The integral equation approach is considered by Lee and 
Bishop [ 10 ], and White [llj. The perturbation method is 
used by Martin [12 ,13] , The frequency response methods are 
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described by Rubin [I4,15]f and Leckie and Pestel [16]* 

It is well knownthat the bending and torsional 
modes get coupled because of the cross-sectional properties. 
In what follows, the techniques generally used for solving 
coupled and uncoupled systems are discussed along with the 
methods to incorporate the other parameters. For the sake 
of convenience this review of the analytical methods are 
subdivided into several sections, in such a manner, that 
in each section a different aspect of the problem is 
discussed. 

2.1 BENDING VIBRATI03!B 

The natural frequencies for cantilever beams of uniform 
crosssection can be determined by classical methods, and are 
given by ^2 

“ fwi; (2.1) 

where El is the flexural rigidity, the mass per unit 
length and L the length of the beam. The values of a^^ 
for the first four modes are 1.875,4.694, 7.855 and 10.996, 
For tapered beams it is more convenient to use 
numerical methods. The tabular method of Myklestad [17 ] 
is, probably, the most widely used. In this method the 
beam mass is Ivunped at a number of discrete stations along 
the beam. A discussion of different methods available for 
f omning the lumped parameter systems is given by Minhinnick 
[18]. Leckie and Lindberg [19] have considered the effect 
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of lumped parameters on bending vibration. 


2.2 fORSIOITAL VIBRAflOJTS 

The classical method for the analysis of torsional 
vibrations of circular sections result in a differential 
equation of motion of the second order. For cantilever 
rods the substitution of boundary conditions of zero twist 
at the fixed end and zero torque at the free end yields 
the expression for torsional frequencies in the form 




niE__v Cg/i 

OT, 


t n — l,3f3,.«.3 

[ 2 . 2 ] 


where C is the torsional rigidity, L the length and I. _ the 

c.g 

polar moment of inertia per unit length of the rod. Sometimes 
this expression is also used to give approximate frequencies 
for non-circular rods by using the value of C given by the 
Saint Venant^s theory of torsion. 

For an aerofoil section the torsional rigidity cannot 
be calculated easily and recourse is generally made to 
experimental and empirical methods. Jawson and Ponter [20] 
have suggested a method for calculating the torsional rigidity 
of any arbitrary cross-section, to an accuracy of one percent, 

Scanlan and Rosenbaum [6] gave Rayleigh’s method, the 
Rayleigh-Ritz method, and matrix iteration procedure for 
calculating torsional frequencies of cantilever blades. Walker 
[21] gave simple formulae for the fundamental torsional mode 
obtained by substituting the original system, by a torsion 
member without mass and having single moment of inertia. 
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Uncoupled torsional vibrations of uniform and nonuniform 
cantilever beams have been considered by Duncan [2], Houbolt 
and Anderson [ 22 ], and ¥ittmeyer [23,24]. Dimcan and Collar 
[3], Pipes [ 25 ] and Thomson [26] gave matrix solutions for 
uncoupled torsional frequencies of cantilever beams. Vet [2?] 
gave results of torsional vibration of rectangular crosssection 
cantilever beams, for general ratios of length to depth and 
width to depth, obtained by the use of the Rayleigh- Ritz 
method . 

2.3 COUPLED BEUDIHU-TORSIOIT VIBRATIONS 

Coupled bending torsion vibrations occur when the 
centre of flexure does not coincide with the centroid of the 
blade crosssection. Garland [28] considered the case of such 
a cantilever blade and used the Rayleigh-Ritz method to 
determine the frequencies and amplitude ratios, Isakson and 
Sisley [29] considered the case of a pretwisted beam vibrating 
in coupled bending.»torsion and determined the frequencies 
using an extent ion of the Holzer-Myklestad method. Rao [30] 
developed a numerical procedure to determine the frequency 
equation of a turbine blade executing bending- tors ion 
vibrations and studied the effects of coordinate distance 
between the centre of flexure and the centroid on the natural 
frequencies of a straight uniform blade with asymmetric 
aerofoil crosssection. Torsion introduces root fixing, i.e., 
at about ( 0.2 X Chord ) of blade length above the root is 

rigid in torsion. 
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2.4 COUPLED BBUDIHrO-BEITDIlTa 7IBRATI0ITS 

If iiie blade has pretwist, uncoupled bending modes 

in flapwise and chordwise directions are mt possible. Both 
the lateral deflections are always coupled and coupled 

bending-bending vibrations occur. 

Dunholter [51] considered the static displacements and 
natural frequencies of pre twisted beams. White [32] employed 
Green's functions to derive the conditions of orthogonality 
for a nonuniform pre twisted blade executing bending-bending 

vibrations, Mendelson and Gendler [33] used a method of 
station fimctions to study the effect of pretwist on lateral 
vibrations of cantilever blades. Diprima and Handelman [34] 
solved the equation of motion of a pretwisted cantilever 
blade by the Rayleigh-Ritz principle, Targoff [35] considered 
the case of a pretwisted rotating beam executing lateral 
vibrations. Carnegie [36] considered the effects of pretwist 
on lateral vibrations of pretwisted blades. Dawson and Carneg- 
ie [37] determined the coupled natural frequencies of a 
pretwisted cantilever blade by a transformation method, Rao 
[58] used Galerkin method to determine the first five matural 
frequencies of a pretwisted tapered cantilever blade, 
Weidenhammer [39] used Variational principle' to compute the 
lowest eigen frequencies of the Qf:)upled system. Using the 
finite difference method, Carnegie and Thomas [40,4l] 
considered pretwisted tapered blading and showed good agreement 
between theoretical and experimental results. 
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The variation in the frequency of vibration of a uniform 
beam due to pretwist depends on the pre twist angle and the 
w id th-to -depth ratio of the cross-section. 

The influence of taper onto an otherwise uniform 
square cross-section beam makes the principle flexural 
rigidities of the beam unequal except at the root. When such 
a tapered beam is pretwisted, coupling occurs between the 
bending motions in the two mutually prependicular planes. 

The degree of coupling is a function of pretwist angle and 
taper, 

2*4.1 TORSIOML VIBRATIONS 

for a pretwisted blade of rectangular section the 
torsion and the bending modes remain uncoupled. The torsional 
stiffness, however, is very much increased due to pretwist, 
resulting in an increase in the torsional frequencies. The 
main cause of the increase in the torsional stiffness over and 
above that given by the Saint Venant's theory, is due to the 
inclination of the blade longitudinal fibres which become 
helical due to the initial twist, 

Chen [ 42 ] considered the effect of pretwist on the 
torsional rigidities. Bogdanoff and Horner [43] gave the 
results for the influence of rotation on the first, three 
natural frequencies in torsion of a unifom cantilever beam with 
different base setting angles and constant pre twist rates. 

Brody and Targoff [44] considered the torsional vibration of a 
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thin pretwisted beam. Carnegie [45] considered an additional 
effect due to torsion on torsional vibration of pretwisted 
cantilever blades. Rao [46] considered the effects of 
fiber bending on the torsional vibration of pretwisted 
cantilever blades. 

2.4.2 COUPLED BEITDIHG-BEMDING-TORSIOR VIBRATIONS 

In a blade of aerofoil cross-section, apart from 
the coupling between the bending modes, produced by pretwist, 
the torsional vibrations are also coupled, Houbolt and Brooks 
[ 47 ] derived the differential equations of motion for 
pretwisted cantilever beams with asymmetrical aerofoil 
cross-section and suggested the Rayleigh-Ritz method for their 
solution. Carnegie, Dawson and Thomas [48] determined the 
natural frequencies and mode shapes for the first five modes 
of coupled bending- bending~torsion vibrations, using a method 
of finite differences. Dawson [49] has determined the natural 
frequencies of cantilever beams executing coupled bending- 
bending-torsion vibrations by Rayleigh-Ritz method. Carnegie 
and Dawson [50] used the method of finite differences and an 
analytical process in which the three coupled equations are 
transformed to an ordinary tenth order differential equation, 
to determine the natural frequencies and mode shapes for a 
straight blade of asymmetrical airfoil cross-section. Rao 
and Carnegie [ 51 ] used Gale'^kin method to determine the first 
five coupled frequencies of a straight uniform blade with 
asymmetric aerofoil cross-section. 
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2.5 EFFECTS OF SHEAR DEFQjRMAJ lOU AMD ROTARY IMIRTIA 

For stubby beams and where higher modes are required, 
the Bernoulli-Buler equation of motion for lateral vibration 
introduces considerable error. Rayleigh introduced the 
concept of rotary inertia, which was later extended by 
Timoshenko [52] who obtained an equation allowing for both 

rotary inertia and shear deformation by a conventional 
Newtonian approach. Carnegie [53] derived the differential 
equation of motion by the use of calculus of variatidns * 

Sutherland and Goodman [54] solved the Timoshenko 
beam equation for the case of simply supported and cantilever 
beams. Chaplin [55] considered the effects of rotary inertia 
and shear deflection while determining the lateral natural 
frequencies of the blades of a 12- stage turbine. .Hurty and 
Rubenstein [56] derived generalized mass and stiffness matrixes 
from the Kinetic and strain energy expressions and presented 
results for a uniform simply supported beam. Lee and Bishop 
[ 57 ] used integral equations for the solution of flexural 
vibrations of a wedge taking into account both the effecis of 
rotary inertia and shear. Carnegie [58] studied the effects 
of shear deflection and rotary inertia for straight and 
pretwisted uniform cantilever beams. Lawson [59] used the 
Rayleigh-Ritz method and showed good agreement between 
his results and those of Reference [54]. Rao [60] used 
the Ritz method to determine the natural frequencies in 
lateral vibrations of tapered cantilever beams allowing for 
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shear and rotary inertia. Krupka and Baumanis [6l] considered 
the bending- bending mode of a rotating tapered pretwisted 
turbomachine blade including the effects of shear and rotary 
inertia. Dawson, Ghosh and Carnegie [62] presented the modal 
curves of pretwisted blading allowing for rotary inertia and 
shear deflection obtained by Runge-Kutta process. Others 
who contributed to the work include Jacobson [63,64], Scanlan 
[ 65 ], Kruszewski [66], Betties [67,68], Mindlin and 
Deresiewicz [69J , Cowper [70], Gaines and Volterra [71] . 

The frequency parameter decreases due to shear 
deflection and rotary inertia as the root width- to- length 
ratio of the beam increases, for any taper* 

2.6 EFFECT OB TABER • 

The mathematical solution for a straight uniform beam 
in simple bending was well known as early as 1736. Meyer [72], 
Ward [ 73 ], Wicholson [74], Wrinch [75], Akimasa [76], Conway 
[ 77 ], Watanabe [78], Eubuo [79], and Wang [ 8 O] gave solutions 
of lateral vibrations of cantilever beams with variable 
cross-section. Southwell [ 8 I] gave a graphical method for 
lateral vibration of nonuniform beams. Matrix methods using 
influence coefficients were given by Duncan and Collar [3], 
Thomson [26], Deckle and Lindberg [82], and lindberg [83]* 
Martin [84] derived correction factors for the effect of taper 
on flexural frequencies by expressing the frequency as *.a 
double power series on two taper parameters of breadth and 
thickness. Taylor [85] solved the basic differential 
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equations of flapwise blade bending by omitting structural 
stiffness terms. Housner and Keigbtley [86] determined the 
natural frequencies and mode shapes of variable cross-section 
cantilever by means of a digital computer, using Prohl method. 
Frequency parameter ratios for the first five modes of tapered 
cantilever beams in lateral vibration are given by Carnegie, 
Dawson and Thomas [48], Rao [87] determined the uncoupled 
bending and torsional frequencies of tapered cantilever beams 
by using the methods of Galerkin. Rissone and Williams [88] 
studied the frequencies of tapered beams over a small range 
of depth ratios. Thomas [89] used the finite difference method 
to determine the bending frequencies of tapered cantilevers. 
Carnegie and Thomas [9Q] gave the natural frequencies of 
long tapered cantilevers using the finite difference method. 

Rao and Carnegie [9l] determined the first three lateral 
frequencies of tapered beams by the use of the Ritz-Galerkin 
process. 

Thomas [89] and Thomas and Carnegie [92] used the finite 
difference method to study the effect of taper on the torsional 
vibration of Cantilever blade. Rao [93] used the Galerkin 
method to consider the effect of thickness taper on the 
torsional vibration of cantilever beams with rectangular 
cross-section. Mabie and Rogers [94] considered doubly 
tapered cantilever beams in transverse vibrations. 

The variation in width taper has little effect on the 
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decrease due to shear deflection and rotary inertia. The 
increase in depth taper increases the effect of shear 
deflection and rotary inertia, thus reducing the frequencies. 

2.7 EPgBGTS OF ROTATION 

The natural frequencies of rotating turbomachinery blades 
are known to ^significantly higher than those of the non-rotating 
blades. The design' speed for a turbomachine is generally 
established by drawing a Campbell diagram [95] for each row 
of the blades. Essentially, this diagram is a plot of the 
variation of the first few. natural frequencies of blades with 
the speed of rotation of the machine, with engine order lines 
superimposed on it. All points of intersection on the diagram 
represent possible resonance if appropriate excitations are 
present. Eor the mechanical design of the blades, therefore, 
the knowledge of the blade natural frequencies, at various 
speeds of rotation, is important. 

Lo and Renbarger [96] derived the differential equation 
of motion of a cantilever blade mounted on a rotating disk 
with stagger and showed that the frequencies of a bar 
vibrating transverse to a plane inclined at an “iwagle to the 
plane of rotation can be found by a simple transformation of 
the frequencies of a bar vibrating perpendicular to the plane 
of rotation. Lo [97] in his analysis simplified the nonlinear 
problem by assuming the blade to bie rigid everywhere except 
at the root and presented the solution in a phase plane. 

Isakson and Eisley [98] considered a similar model to 
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determine the natural frequencies in coupled bending and 
torsion of twisted rotating blades. Boyce, Diprima and 
Handleman [99] used Rayleigh-Eitz and Southwell methods to 
determine the upper and lower bounds of natural freq,ueiicies 
of a turbine blade vibrating perpendicular to the plane of 
rotation. Yntema [lOO] making use of the bending mode 
of a nonrotating beam determined the first three modes of 
a rotating beam by the use of a Rayleigh energy approach.. 

Boyce [lOl] used Rayleigh-Rttz and Southwell methods to 
detennine the upper and lower bounds of natural frequencies 

and showed that the frequencies depend almost linearly on 
the hub radius for various rotational speeds. Carnegie [L02] 
derived a theoretical expressions for the work done due io 
centrifugal effects for small vibrations of rotating cantilever 
beams and established an equation for the fundamental 
frequency of vibration by the use of Rayleigh's method. 
Schilhansl [l03] investigated the stiffening effect of the 
centrifugal forces on the first mode frequency of a rotating 
cantilever blade by the use of successive approximation, 
Carnegie, Stirling and Fleming [104] solved the 
differential equation of motion by the finite difference 
method to study the effects of rotation and stagger angle of 
the blade. Kissel[l05] has given the modified Southwell 
factor for the effect of rotation, Rao [106] used the 
Galerkin method to study the effects of rotational speed, disk 
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radius and stagger angle of the blade and obtained a simple 
relation to determine the frequency of lateral vibrations. 

Rao [107] developed a numerical procedure to determine a 
polynomial frequency equation of a rotating cantilever beam. 
Rao and Carnegie [108,109] used the Ritz energy process and an 
averaging, procedure to solve the differential equation derived 
by Carnegie [110] and obtained the nonlinear response of a 
cantilever blade. Ansari [ill] applied Perturbation and Ritz 
energy methods to the problem of non-linear flexural vibration 
of a twisted, rotating cantilever blade of ssrmmetrical cross- 
section vibrating in the plane of rotation. Pnueli [ll2] 
compared natural bending frequency to rotational frequency 
of rotating cantilever beam. Kundu [113] considered 
transverse vibration of an isotropic rotating, variable cross- 
section beam. Wauer [II4] determined the eigenvalues of 
rotating turbine blade by means of Ritz method. Trompette 
and Lalanne [115] did the theoretical and experimental 
vibration analysis of rotating turbine blades, including 
large strains, root's influence, temperature effects using 
the finite element method. Whereas, Henry and lalanne [II6] 
performed the vibration analysis of rotatiJig compressor 
blades. Ansari [ll7] .presented an analysis for evaluation 
of non-linear modes of vibration of a pretwisted non-uniform 
cantilever blade of uns 3 nnmetrical cross-section mounted on the 
periphery of a rotating disk, including shear deformation, 
rotary inertia .cariolis effects, using the discrete system 



19 


formulation. Others who contributed to the work include 
Liebers [118,119], love and Silberstein [120], Billington 
[121], Sutherland [l22], Kobori [l23], Plunkett [l24 ] , liner 
[125] and Griansante [126], 

2.8 EFFECTS OF ROOT FLEXIBILITY 

In most of the theoretical investigations, the blade is 
assumed to be completely fixed at the root, in practice, 
however, the blade is rather loosely mounted on the rotor 
disk* The vibrations of the blade, therefore, continues to 
some extent into the mounting. This has an effect similar 
to that of increased blade length and results in a lower 
frequency than the one calculated on the fixed end assumption. 
The root flexibility, generally, makes the root a non- 
deflecting end, but not a non-rotating end. As a result, the 

slope at the root is not zero but is proportional to the 
moment acting at the root. The root flexibility varies with 
the speed of rotation, since the centrifugal forces tend to 
tighten the blade in the disk and increase the fixidity. 
Moreover, various types of root fixings are used in practice 
and, therefore, a general analytical treatment of this effect 
is rather difficult. Some investigations, however, have been 
carried out and are discussed here. 

Niordson [l27] has investigated the vibrations of 
turbine blades with loose hinge support in centrifugal force 
field. Hirsch [l28] considered the effect of non— rigid support. 
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Traupel [129] has suggested an empirical method of taking 
into account the effect of root flexibility on the fundamental 
natural frequency in bending. He proposed the relation, 

where is the frequency calculated on the fixed end 
assumption andi^jj^ is the actual frequency. He has plotted 
the value of against the ratio (^ ) of the blade length 
to the radius of gyration of the root section, which is 
considered as a measure for the flexibility of the root. Baur 
[130] has suggested that the value of given by Traupel [129] 
may also be used for the second flexural mode. For torsional 
vibrations, Baur [130] has suggested a similar relation, but 
has pointed out that the variation of the correction factor 
be plotted against the non-dimensional blade flexibility Z /G, 

Jr 

where I and C are the polar moment of inertia and the 
^ at 

torsional stiffness/the root, respectively. Horway [l3l] 
considered hinged ends and Bogdanhoff [l32] considered the effect 
of pinned ends on the natural frequencies of rotating blades. 

The problem of pin— fixed blade is considered in more 
detail by Goatham and Smailes [133] . The blade is considered 
mounted on the disk through a loosely fitting pin, with a 
clearance. The pin tightens when the blade rotates, due to 
the centrifugal forces on the blade. Shaw [l34] proposed 
an improved blade root design for the axial flow compressor 
and turbines. 
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2.9 EFFECTS OF DISK ELASTICITY 

Jet engine raters are highly stressed. Considerable 
care must be taken for designing such systems to ensure that 
the parts do not fail during its life time. 

If the disk on which the blades are mounted is 
considered absolutely rigid, each blade can be treated as 
a separate ixnit for vibration analysis. However, the disk 
has some elasticity, and hence the blades and the disk ought 
to be considered as an assembly* Energy from a blade can be 
transferred to the adjacent blades through the disk, which 
modifies the vibration characteristics of the blades. 

Ellington and McCallion [l35] have investigated the 
effect of elastic coupling through the disk rim, on the 
frequencies of vibration of the blades, by considering a 
simplified model. In this model, the blades are replaced 
by uniform beams fixed to the rim at their roots and vibrating 
in a plane parallel to the plane of the disk. The rim is 
considered as a uniform elastic ring, permitting no radial 
displacement at the roots of the blades. The mass of the 
rim is either neglected or considered as being concentrated 
at the root of the blades. The analysis has shown that if 
there are H identical blades mounted on an elastic rim, 
the system has N fundamental frequencies, H first overtones 
etc. However, if stiffness of the rim is large compared to 
that of the blade, as is generally the case, these H 
frequencies are very near to the corresponding frequency of 
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the individual blade. Thus instead of individual natural 
frequencies, there are hands of frequencies — width of hand 
decreasing with the increasing stiffness of the disk. This 
explains to some extent the scatter in the natural frequencies 
usually observed in vibration tests. 

The coupling effects due to the blade ring, for the cases 
when either all blades have identical frequencies or alternate 
blades have identical frequencies, have been investigated by 
Sohngen [l36]. Pillipov [l37] has considered the problem 
of the tangential vibration of the blades together with the 
torsional vibrations of the disk, when all the rotating blades 
vibrate in phase. A simplified model for a row of blades 
mounted on a flexible disk is suggested by Wagner [138]. By 
using a model of blade-disk assembly more or less similar 
to one suggested by Wagner [138], Dye [l39] has studied the 
effects on the stress level of the blades due to different 
factors. The factors taken are ; 

(a) the blades of a particular stage having slightly varying 
natural frequency and, (b) the transmission of forces from a 
blade to its neighbour through the flexible disk. 

Shen [l40] developed a mathematical formulation for the 
analysis of transient and steady state flexible rotor 
dynamics, including the general non-axisymmetric and non- 
synchronous rotor motion which may result from the included . 
in -phase and out of phase stiffness and damping functions at 
all bearing and rotor stations. Other rotor dynamic parameter 
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considered are the rotor masses and mass moment of inertia 
and their eccentricities and misalignments, 'The effects of 
rotor drive and disipative torque and the interaction between 
tbrsional and transverse motion are also included, 

Lemke and T rumple r [l 4 l] presented an analytical 
investigation of the effect of coupling characteristics on the 
response of axially coupled sets of turborotors. The approach 
used is transfer matrix technique which utilizes a generalized 
machine element matrix and coupling fixity matrix. Disk 
inertia, gyroscopic effects, and shaft as well as coupling 
orthotropy are included in this study. 

DopKin [142] developed a general analysis to predict 
the natural frequencies of vibration of a rotating shaft 
on linear, un-fdamped bearings j including flexible disks. The 
relative merits of several possible analytical methods are 
examined, and the transfer matrix technique is chosen to 
solve the system. The shaft and disks are described by a number 
of constant thickness segments and lumped masses, A fortran 
IV computer program is presented which solves the system 
equation by successive trials. Results of the analysis are 
found to agree very well with the results obtained by others 
for either flexible disks or flexible rotors with stiff disks, 

A limited parametric study of six different rotor 
geometries with a range of disk flexibilities is presented. 

Some of the conclusions resulting from this study are : 
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(1) Disk flexibility always reduces the rotor natural 
frequency, sometimes to a few percent of the rigid disk 
value . 

(2) When disk frequency is greater than four times rotor 
frequency, disk flexibility has a negligible effect. 

(5) When the disk and rotor have the same frequency, the 
combined frequency is in the neighbourhood of 2/3rd as 
great, 

(4) When rotating speed is greater than l/2 of the natural 
frequency, a lower limit on natural frequency can be obtained 
with a rigid disk analysis by setting the disk polar and 
diametral moments of inertia equal to zero. 

2.10 EFFECTS OP PfiETWIST 

If the blade has pretwist, uncoupled bending modes 
in flapwise and chordwise directions are not possible. 

Both the lateral deflections are always coupled and coupled 
bending -bending vibration occur. 

Young [143] presented the effect of twist on blade 
vibrations. Thomas and Dokumaci [l44] used Finite element 
method for pretwisted blade vibration. From the results 
obtained by various investigators, certain general conclusions 
can be drawn regarding the effect of pre twist on the frequencies 
of cantilever beams of rectangular cross-section. The 
conclusions are valid for the magnitudes of pretwist usually 
encountered in ttirbomachinery blading : 
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(a) The fundamental frequency is affected -vB^ry little by 
the pretwist. 

(b) If a chordwise frequency and a flapwise frequency are 
very close to eachother, the effect of pretwist is to separate 
them more, 

(c) If the depth to width ratio is very large (say more than 
20) the second and the third bending frequencies decrease 
with pretwist and 

(d) The effect of twist on the natural frequencies is less 
in tapered beams than in the beams of uniform cross-section. 

2.11 EFFECTS OF TEMPERATURE 

The first effect of the high temperatures is a deteriorat 
ion of the materials of blading. Hoff [145] has collected a 
great deal of information on the variation with temperature 
of the yield stress, the ultimate stress, and Young’s modulus 
of elasticity. 

An equally important second effect, which can often be 
observed at temperatures that do not impair the allowable 
stresses, is the development of thermal stresses. G-ache and 
Ronzin [146 J have shown that the magnitude of vibratory 
stresses is influenced not only by the magnitude of the 
quantitative magnitude of the circumferencial non-uniformity 
of turbine inlet temperature field but also by nature of the 
non-imif-ormity. G-orelkin and Bogov [l47j considered a solution 
of the thermo elastic problem for gas turbine rotor blades 
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and presented the results of theoretical and experimental 
investigations of thermoelastic stresses in these blades* 

A third effect of the high temperatures is that the 
material begins to flow slowly under comparitively small 
stresses j this flow is designated as creep. In a complex 
structure the stress distribution is governed by the equations 
of- equilibrium and the conditions of the compatibility of 
the deformations. ¥hen the latter are established with the 
aid of Hooke's linear law of elasticity, the ensuing stress 
distribution is entirely different from the one that prevails 
when the connection between stress and strain is griven by a 
highly non-linear creep law.. 

Moreover the creep deformations are dependent on time. 
Depending on the material, stress, and the temperature, 
the time elapsing between the load application and fracture 
the so called critical time, may be a few years, a few days, 
or even a few seconds. 

Turbine blades are subjected to a more severe and 
complex combination of stress and temperature than any other 
jet engine component. The blades in current engines are subje 
cted to very high centrifugal stresses at temperatures of the 
order of 1200°C. Also they are subjected to hot corrosive 
atmosphere and to gas impulses that may cause the blades to 
vibrate. Rapid • beating and cooling induces thermal stresses 
in them. Hence turbine blade failures are very common and 
occur in a variety of ways. 
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In addition to the centrifugal load, the gas forces 
impose a bending load on the blade airfoil. The designer 
partly compensates for the gas load by tilting the airfoil 
slightly do¥ 2 istream so that centrifugal force will induce 
opposite bending. The bending force can be completely 
cancelled at only one plane along the airfoil span. Further 
since the gas bending load reduces with altitude and engine 
operates at essentially constant speed, the gas bending load 
can be cancelled for only one altitude. If the blade is 
tilted to compensate for the sea level gas load, it will be 
overtilted for the gas load at altitude. 

2.11.1 STRESS-RUPTURE 

In the range of temperatures in which turbine blade 
operate, there is, for each blade material, a finite time 
before fracture for each combination of stress and temperature. 
In general, the life predicted from stress -rupture properties 
is the longest that the blade could be expected to run at full 
engine power, since. the only stress considered is that 
induced by a centrifugal force. Other environmental factors, 
such as vibratory stress, corrosion, thermal stress and impact, 
will reduce the life belovj this value* 

Incipient stress -rupture failure appears as irregular 
intergranular cracking in a narrow zone of airfoil span known 
as. ’ '.Critical zone Usually cracks are not confined to 

the leading or trailing edge, but occtir at random across the 
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chord. Blades that have completely fractured have many cracks 
on the airfoil surface adjacent to the fracture. Since the 

fracture is intergranular, the fracture surface will generally 
he rough, 

Linask [l48] conducted an analjijical study using 
Fracture maohanios principle to model turbine airfoil cracking. 
He found that crack initiation can be related to calculated 
residual strains in the airfoil coating and that coating 
properties are an important consideration in determining crack 
location and orientation. Ihe coating crack subsequently 

propagates into base material according to basic fracture 
machanics laws. 

Detection of stress -rupture cracks during engine 

inspection cannot be counted upon as method of avoiding blade 

fracture, because complete fracture generally follows cracking 
very shortly, 

2.11.2 nRW-RP 

At any condition of stress and temperature on the 
stress rupture curve, the test specimen elongates with time 
in a manner described by the ideal curve of Fig. (2. 11,2 ), 

When load is applied, the specimen elongates elastically 
then plastically at a decreasing rate (first-stage creep) 
until the rate becomes approximately constant (second stage). 

Finally the rate begins to increase (third stage) until 
fracture occurs. 
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Unit elongation or strain, is not uniform along ■ 
a blade length, because of the nonuniform stress and 
temperature conditiondsTfee zone of maximum creep corresponds 
to zone of minimum stress - rupture life. 

The designer is interested in knowing the allowable 
centrifugal stress and temperature because they affect both 
the time to rupture and the time in which elongation of the 
blade will deplete the operating clearance and cause the 
blade tip to rub shroud band. Because the strain in a blade 
is highly localized, total elongation is usually small t stress- 
rupture fracture usually occurs before the allowable total 
elongation is exceeded. Thus, stress -rupture usually becomes 
the more important design criterion. 

Creep appears as blade stretch, blade- shroud rubbing, 
or blade necking (locMl reduction in area). Necking cannot 
be relied upon to indicate incipient blade fracture because 
the time between its first appearance and fracture is usually 
too short and unpredictable. 

2.11.3 THEBMAL STRESS FATIGUE 

large temperature gradients occur in turbine nozzle 
vanes and blades during transient conditions such as start| 
up , acceleration, and shutdown. Since the leading and 
trailing edges of turbine blades heat and cool much more 
rapidly than the thick midchord region, chordwise temperature 
differences occur in the blades whenever the engine is operated 
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i/hrough. "brans ienijs. For example, in a normal S'tar't on one 
engine, a difference of 500°C was measured between the 
leading edge and the centre of a turbine blade about 8 seconds 
after ignition. 

Temperature differences induce stresses in the blades 
because of differential thermal expansion. The heavier body 
section of the blade restrains the expansion of edges, and 
thus high compression stresses (during heating) or tensile 
stresses ( during cooling) are induced in the edges. With 
the result that plastic flow may occur in these areas. When 
the body eventually attains the equilibrium temperature, the 
direction of the thermal stress in edges which have been 
plastically deformed is reversed. This phenomenon is repeated 
during successive cycles of engine operation and has caused 
warping and cracking in edges of nozzle vanes and blades. 

2 .12 FATIGUE 

Periodic forces resulting principally from disturbances 
in the gas flow through the engine cause vibration of such 
items as turbine blades, compressor vanes and blades. When 
the frequency of the vibratory force is in resonance with 
a natural frequency of a component and the damping losses 
are small, high stresses can result which may cause early 
failure of the component. 

2.12.1 COMPRESSOR BLADES FATTGTTE 

One of the most vulnerable components in a gas-turbine 



31 


unit is the axial -flow compressor. Vibration fatigue account 
for many compressor losses. 

Fatigue usioally occurs in the rotor blades, but is 
occasionally experienced ani the stator blades or rotor 
disks at the blade recess because of vibratory loads from 
the blades. The blade vibrations are generally induced by 
pulsations in the airflow. In most cases of blade failure, 
the fractures occur in the airfoil section near the base and 
show progressive damage typical of pure fatigue. The origin 
of the failure is usually on the convex surface at the maximum 
camber point, where the bending stress is maximum during 
vibration. Fatigue failures may also originate in the trailing 
edge of the airfoil. In other cases, the failure is in the 
fastening. Frequently fatigue failures originates at random 
locations on the airfoil where stresses are not ordinarily 
high but are magnified by stress concentrations arising from 
toolmarks or nicks made by foreign objects. Six types of 
fastenings have been used for compressors. In the fir - tree 
types, high stress concentrations occur in the fillet between 
serrations, particularly if the radii are small. The bending 
stresses due to vibration are superimposed on the steady 
centrifugal and gas-bending stresses. All stresses are 
magnified by the stress concentration caused by the abrupt 
changes in cross-section. 

Fatigue results from alternating stresses produced by 

. listed in order of importance, factors that 


blade vibrations 
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may induce these vibrations are : 

(1) Rotating stall, (2) aeroelastic coupling. 

2.12.1.1 ROT Al IRQ STALL 

Rotating stall is one of the principal causes of seVere 
vibrations leading to fatigue failures . This phenomenon 
occurs mainly during acceleration of the engine through 
the speed range from 50 to 70 percent maximum speed. The jet 
engine is seldom subjected to this speed range for long 
periods of operation, as the range is traversed only upon 
acceleration and deceleration. Neverthless, rotating stall 
constitutes the primary source of faitgue failure. Reference 
[ 149 ] gives means of alleviating rotating stall frequencies. 
These include (l) increasing the damping of the blades, (2) 
preventing the operation of any compressor blade row in a stall 
ed attitude, and (3) disrupting the periodicity of the stall 
patterns so that a resonant condition between the blade 
frequency and stall frequency cannot exist. 

2.12.1.2 AEROELASTIC 00UPLIR& 

The interaction of aerodynamic forces and the elastic 
behaviour of the blading can result in an Instability which 
produces blade vibrations. Flutter is an instability- 
self maintained by the continual absorption of energy from 
the airstream. In compressors, the problem may be aggravated 
by high oscillatory aerodynamic loading, proximity to stall, 
pressure rise across the stage, and certain cascading effects. 
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Excitations occuring in the vicinity of stall is known 
as Stall flutter, which is entirely different from rotating 
stall. Rotating stall is a periodic nonuniformity in the air 
flow that only excites blade vibration when the impulse 
frequency approaches the blade natural frequency [ 150 ] * 

Stall flutter is self - induced vibrations of the blade in a 
xmiform flow field due to aerodynamic instabilities arising 
from high angles of attack. 

Stalling flutter has been shown to occur on compressor 
cascades in bench tests [l5l] . Experimental investigations 
[ 152 ] carried out on compressors under operating conditions 
have shown that in some cases the vibratory stresses increase 
greatly with angle of attack. An analysis was made of the 
effect of centrifugal forces on the flutter of a uniform 
cantilever beam at subsonic speeds with application to 
compressor and turbine blades by Mendelson [153]* This is 
dealt in more detail in section [ 2 , 13 ]. 

2.12.2 TURBIEE BLADES EATIGUE 

The most important source of excitation of turbine 
blade vibrations is the impulse given to the blades when 
they pass through the wakes of the nozzle vanes. The flow 
from individual combustion chambeis likewise introduces 
irregularities into the gas flow impinging on the blades 
and may also excite vibrations. The frequency of vibratory 
impulse imparted to the blade is thus a function of the 
number of nozzle partitions, or cumbust ion chambers, and of 
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engine speed. The blade can vibrate in simple bending, 
in torsion, or in complex combination. 

When the vibrating stress is superimposed on the 
centrifugal stress the appearance of the fractures is 
dependent upon the magnitude of the vibratory stress relative 
to mean tensile stress. This is shown in reference [l54]. 

A combined creep and fatigue problem may arise in aircraft 
gas turbines, where the turbine blades are subjected to centrif- 
ugal stresses together with fairly high frequency (1000 Hjs) 
bending fatigue stresses caused by the hot gases impinging 
on the blades j temperature may be above 1000°C « The rotor 
disks are subjected to repitetion of tensile stress for 
several hours, together with p'eak stresses for short times. 

High local stress concentrations are set up at cooling holes 
and fir - tree blade root fixings, and the material must be 
capable of with standing these for approximately 10000 fli^t 
cycles. It is stating the obvious when it is noted that 
the failure of blades or disks may have disastrous effects 
on the safety o,f the aircraft. Ellison [l55 ] lias discussed 
modes of failures and various methods of testing where 
there is . ihteraction between creep and fatigue# 

Majumdar [156] has analysed all the major components 
of the gas turbine engine for low cycle fatigue and /or creep, 

A linear cumulative damage law is assumed, since the 'interaction 
between creep and fatigue is still an unsolved problem. 
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2.13 AERODYNiHIC EXCITATIONS 

Tlie discussion so far was confined to the methods 
for the calculation of natural frequencies and mode shapes, 
however, from strength consideration point of view, knowledge 
of the amplitude of vibration is also essential. The aerodynamic 
exciting forces, the methods of calculating deflection 
amplitude, stress levels, etc. are considered. 

The aerodynamics of compressors, especially axial 
flow, being complex, the analytical treatment of this field 
is somewhat limited. The practical aspects of design and 
vibration analysis have been reviewed in several papers. The 
more important of these reviews are those by Blackwell [l57], 
Pearson [158], Armstrong and Stevenson [159], Ref.[l54] and 
Ref. [149] etc. 

2.13.1 FLUTTER 

The approach to flutter prediction of fans and compressors 
in gas turbine engines is different from the approach to 
flutter prediction of aircraft wings and other external 
aerodynamic systems. These differences arise from the greater 
disparity between the density of the air and the density of 
the structure, and from the more complex boundary to the air 
formed by adjacent blades and walls in the gas turbine. 

In the past [Ref. 149,154,157,158,159], the diffictilties 
of predicting the proper tlnsteady aerod 3 mamic forces for 
cascades of blades found in turbomechines have necessitated 
the use of empirical approaches to flutter prediction. These 
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approaches tiled to distinguish many different types of 
flutter and to correlate parameters measuring the onset 
of flutter with other structural or aerodynamic parameters 
[l60-l63]» Thus the onset of flutter was usually characterized 
by the relative velocity at some specified spanwise location 
on the blade, divided by the product of vibrational frequency 
and chord of the blade at blade at that spanwise position. 

This was correlated against a number of parameters like the 
ratio of torsional to translational deflections, the relative 
Mach number,, or the nondimensionalized incidence angle. The 
parameteric plots served to correlate experimental data for 
different types of flutter [164-166], the onset of flutter being 
defined by an empirical boundary. 

The empirical correlations have become inadequate as 
recent advances in technology (higher tip speeds in turbines, 
avoidance of heavy shrouded construction, etc. ) required 
extrapolations of past experience to new areas. The sharp 
dividing lines between stable operation and flutter disappeared. 
Attempts to improve the correlating parameters were unsuccessful, 
due to scarcity of new data and the necessity to extrapolate 
outside past experience. A need to develop a unified flutter 
prediction based on basic aerodynamics and vibration analysis, 
and capable of application to all new design, became very real. 

A uniform approach to flutter prediction has been 
developed by Mikolajczak et al. [167]. The aeromechanical 
stability of the blade -disk system is expressed in terms of a 
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stability parameter which, measures the amount of unsteady 
work done by the air on the systam, when the system is 
vibrating in one of its natural modes. In neutrally stable 
systems, the unsteady work done by the air on the blades 
will balance the work dissipated by friction and by material 
damping. An accurate prediction of the vibrational deflections 
and of the unsteady aerodynamic forces is req.uired at eveiy 
spanwise location on each blade, so that the work done by the 
unsteady aerodsmamic forces may be calculated. Recent progresss 
is described in the prediction of unsteady aerodynamic forces 
and the determination of mode shapes. The stability model is 
applied to the prediction of supersonic flutter, chordwise 
bending flutter, and stall flutter. 

Plutter accurs when the unsteady aerodynamic forces and 
moments created by periodic blade vibrations do positive 
aerodynamic work on the blade during each vibration cycle 
and the mechanical damping is insufficient to dissipate this 
work input. 

If the rotor system (blade, disk and shroud ) is 
assumed to be a simple linear system, then the aerodynamic 
damping for system [l67] can be written as 

6 = - nW /4 S’ 

aero 

where n is the number of blades in the rotor, and EE is the 
average vibrational kinetic energy of the rotor system. 

Flutter is obtained when the total damping for any 
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■vibrational mode is less than zero, i.e,, 

64. . T = 6 +6 , '*^1 o 

total aero mech ^ 

The mode "with the lowest damping is assumed to be ihe 
least stable. For flutter prediction, it is usually sufficient 
to examine only the lower frequency modes, namely the first 
three modes, since higher modes tend to have a hi^ positive 
damping. Of course, it is necessary to examine the higher 
modes, so called plate modes, using chordwise flutter 
predictions. 

Boundaries for the most common types of flutter are shown 
schematically on a compressor map in Fig. [2.13.1]. The 
distinguishing features of each type of flutter prediction 
are sximmarized as follows [I68], 

(1) Supersonic unstalled flutter ; This type of flutter 
occurs when the outer span of the blade is operating supersoni- 
cally and can occur at the design point or at either higher 

or lower pressure ratios. Its occurence imposes a limit on high 
speed operation. This type of flutter usually occurs in 
vibration modes which are predominantly torsional [165]* 

(2) Subsonic stalled flutter ; This type of compressor flutter 
is encountered near the surge line . In a high-speed fan stage, 
subsonic stalled flutter may occur at part speed ; in a low - 
or high - pressure compressor, subsonic stalled flutter may 
occur at or near design speed. It is suspected that during 
stalled flutter, separated flow may exist on the compressor 
blades. 
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(3) Choke ( ne gat ive incidence) f lut t e r : This type of 
compressor flutter may occur when a fan is operating either 
suhsonically or transonically. The stresses will increase 
as pressure ratio is lowered* During choke flutter, the rotor 
is operating at near - choke conditions and the flow is 
transonic over most of the blade chord. 

(4) Supersonic stalled flutter : This type of compressor 
flutter occurs when the outer portion of the blade is operating 
supersonically, and the stage is operating near the surge 
line. The flow has strong shocks in the blade passages. 

Stalled flutter and supersonic unstalled flutter are 
the two most common types of flutter. 

Chordwise Blade Flutter 

For very high speed fans for which composite materials 
are currently being developed, supersonic unstalled flutter 
can occur in a vibratory mode which involves a large 
vibrational deformation of the blade cinamber line. Such modes 
are present and important near blade tips when the rotor 
aspect ratio is low and the blades are thin. It can be seen 
that the susceptibility to flutter increases with Mach nmnber. 

Stalled Flutter 

Stalled flutter usually occurs near the compressor 
surge line and is identified by an increase in flutter 
stress as pressure ratio is increased at constant speed. 

In some instances the stall flutter boundary has been 
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observed to occur near maximum compressor efficiency indicating 
that stalling is not an essential condition for ''stalled'* 
flutter, although it could he the most severe condition. 

This indicates that the influence of flow incidence on the 
unsteady aerodynamic forces will have to he considered in the 
prediction of the flutter hoimdary. 

2.13.2 WAZE EXCITATION 

Rotating blades provides some of the most challenging 

problems of fluid mechanics. Plow within a turbomachine is 
complex, three dimensional, rotational, asymmetric, compressible 

viscous and unsteady. 

Wake excitation phenomenon gives rise to an unsteady 
pressure distribution and lift force on each blade is a 
source of mechanical vibration as well as acoustic- 
radiation. The problem is essentially one of unsteady airfoil 
theory, i.e., an airfoil travelling through a sinusoidal gust 
(velocity perturbation). Naumann and Yeh [169] develops 
the unsteady pressure and lift for an airfoil that has both 
camber and angle of attack, under longitudinal as well as 
transverse velocity perturbations, ignoring the interference 
effects due to adjoining blades and hence is valid for cascade 
of low solidity. 

Physical mechanisms (e.g. self - excited vibrations of 
rotor blades and stator vanes, shaft and disk vibrations etc. ) 
are discussed [170] as possible aeromechanical responses to 



41 


distorted flow, Tlie role of various forms of damping and 
the use of composite materials is also described. 

Denforth [l7l] has described 8 mechanisms of distortion 
induced vibrations. A distortion index for blade vibrations 
is defined both as an initial vibration design alert and for 
automatically assuring integration of design for inlet 
aerodynamics with that for blading. These 8 mechanisms are ; 
(1) stall, (2) resonance , (3) random vibration generated 
by excitations associated with cascade separation, (4) its 
combination with resonance maxmized at a critical operating 
temperature, resonance indirectly generated by distortion 
accentuation of these mismatched stator cascades, (5) blades 
instability, or flutter, induced by distortion in stages stable 
without distortion, (6) inlet gust excitation, (7) upstream 
tiirbulence, and (8) oscillatory inlet flows. Discussion is 
limited to blade vibration in the prestall operating regime; 
of fans and compressors ; for stall, referred to fans and 
compressors as a whole rather than local ’’cascade stall’’, 
is a generally unacceptable operating conditions. Blade 
response in fan or compressor stall is a specialised subject 
in itself, 

2,14 ASPECT RATIO 

Much of the early work carried out in studying the 
vibration characteristics of turbine blading was based on 
the assumptions of large aspect ratio, Carnegie [36], Esakson 
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and Eisley [29], and Houbolt and Brooks [47] » among others, 
used simple beam theory and various solution techniq.ues to study 
the vibratory behaviour of twisted cantile'^®^ beams as an 
analog for twisted propellers and turbine hLahing. While 
these works are complete in their own rigl^’*"» their techniques 
can not adequately handle the newer lightweight, lo'W aspect 
ratio turbine blading where the blades are more likely to behave 
as plates or shells rather than beams. The solution to this 
dilemma arrived with the advent of finite element analysis 
methods. Eow, a blade having a complex geometry can be 
accurately modeled as a composition of many elastic elements 
having well-defined elastic properties. Anderson, Irons, and 
Zienkiewicz [l72] aptly demonstrated the utility of finite 
element analysis in studying the vibration and stability 
behaviour of thin plates for various geome't^^®® boundary 

conditions. Rawtani and Dokainish [l75»l''^'^^ have very 
ively used finite element techniques to study the vibration 
behaviour of twisted cantilever plates. They analyzed the 
effect of varying degrees of pretwist on the natural o . 
and frequencies of cantilever plates [175] • They al o 

carried out a vibration analysis of a flat cantilever p a 
fixed to a rotating disk [174]. The effects of plate inclination 
disk radius, and angular velocity on the plntes nat 
frequencies and mode shapes were studied. 

Dokainish [l75] has extended the transfer matrix 
technique generally known as the ' 'Holzer-^yklest ad” method. 



to reduce the size of matrices in the finite element method 
used in the vibration analysis of plates and shells* As 
such the transfer matrix has been applied eahlier by Leckie 
[ 176 ] to plate vibration problem. 

Macbein [177] have presented a combined numerical -* 
experimental study of the effects of yarying tip twist and 
increasing centrifugal loading on the resonant characteristics 
of cantilevered plates. The finite element computer program 
HASTRAN, is used to compute the natural frequencies, mode shape 
and normalized shear stress distribution. 

2.15 EFFECTS OP SHROUB 

Turbomachinery blades are often joined together, either 
at their tips or at some intermediate location, to reduce 
the resonant amplitude. A single shroud ring may join all 
the blades in a row or else a packet of blades may be joined 
together. The continmous shroud is broken between groups 
of blades with a gap of order of 0.5 mm to allow for its 
expansion. The shrouding significantly alters the nat\iral 
frequencies and mode shapes and hence the packet of shrouded 
blades should be analysed as a unit for the tangential, axial 
or torsional vibrations. The shroud also causes a coupling 
between the axial flexural vibrations and the torsional 
vibrations. 

The effect of shrouding on the tangential vibrations 
(vibrations in the place of the disk) of a packet of blades 
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is investigated ty Smith [178] using several simplifying 
assumptions. Prohl [179] considered not only the tangential 
vibrations but also the axial and the torsional vibrations, 

A method of finding the natural frequencies of a group of 
shrouded blades, using the perturbation procedure is given 
by Tuncel et al [180] . The method is applicable for weakly 
cQuple<( blad e and shroud subsystems. Pugino [l8l] has 
considered the effect of shrouding by taking the vibration 
form of the blades as a linear combination of fixed-free, 
f ixed-supported and fixed-fixed modes for the tangential 
vibrations. 

The problem of two adjacent blades, connected at 
their tips by a shroud to form a continuous frame, has been 
analysed by Singh and Nandeoswaraiya [l82] . The fundamental 
frequency is obtained for the tangential and the axial 
vibrations. For rotating blades, a Southwell type correction 
factor is? also derived. 

An investigation of the vibrations of a coupled disk- 
blade-shroud assembly is considered by Stargardter [183], by 
experiments on flexible silicon rubber models. It iMt 
observes that the disk vibrates as a plate, mainly with 
modal diameters, and the shroud acts gs ring. Bhide [l84] 
and Bajaj [185] have also studied the vibrations of packetted 
turbine blades. 

Shroud design is effective against untwist [186], Not 
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only turbine blades, but almost all of tbe large fan blades 
have part-span shrouds to get better vibration characteristics 
and these shrouds stand against untwist deformation by 
constraining themselves eachother* 

2.16 LACINQ WIRES 

One of the reasons for using shroud is to hold blades 
in the correct position. This is also achieved by lacing 
or Lashing Wires. Shrouding introduces high centrifugal forces 
where as Lacing achieves the purpose of alignment without 
appreciable addition to centrifugal forces. 

Formerly lashing wires were strung throu^ the bossed 
holes, machined in the blades. The present practice is to 
weld short pieces of lashing wires to each side of the blade. 
When the blades are machine-finished and stresses are relieved, 
these lashing wires touch each other and are joined either 
by a weld or by a welded sleeve. The main disadvantage of 
lashing wires is that it disturbs the flow pattern and 
contributes to vibration tendencies. To avoid these, airfoil 
section is used for the lashing wires. Lashing wires could 
be provided at two or three places along the length of the 
blade in order to give it sufficient rigidity. 

2.17 LARGE AMPLITUDE 

All most all the papers discussed earlier assumes that 
amplitudes are small. In dealing with large amplitude 
vibration of deformable bodies such as beams, plates and shells 
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■fche relatiionsliip between extensional and shear strains 
on the one hand and the displacement components on the other 
are assixmed to be non-linear. Thus we have geometric non- 
linearity and as a consequence the governing differential 
equation will turn out to be non-linear. This is true inspite 
of the fact that the relations between curvatures and 
displacement components are assumed to be linear, (This is 
also called as Structural non-linearity — due to large amplitr- 
ude causing stiffening of the structure thus raising the 
frequency of response ). 

This need not be the case if the non-linearity is due 
to material ( e.g. Ramberg- Osgood material). 

Another way of looking at this is to consider the total 
strain energy of the deformable elastic body which is the 
STim of the extensional or stretching energy and bending 
energy. The extensional energy becomes non-linear (involving 
higher order . non-linear terms than quadratic) in the normal 
displacement component, whereas the bending energy remains 
quadratic in the displacement components. 

Srinivasan [l87] has considered non-linear vibration 
of beams and plates, Saythyamorthy and Pandalai [l88] has 
obtained governing equations for large amplitude flexural 
vibrations of plates and shells, large amplitude vibration of 
deformable bodies are considered in Ref, [189] • 

2.18 DAMPm 

The damping in a system determines the amplitude, and 
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conseciuently the vibratory stress, for a given exciting 
force. A material may have very high fatigue strength, 
but if the damping is very lew, the material may have a much 
shorter life than one with moderate fatigue strength and 
high damping for the same exciting force. Material damping 
is the most reliable of all the sources available. Root 
damping may vanish because of high centrifugal force or 
malfunctions of an external damper. Aerodynamic damping 
may vanish because of unfavourable coupling betvreen the 
aerodynamic and elastic forces. Some materials display 
increased internal damping as the vibratory stress increases 
Ref. [154]. 

Various types of damping can be employed ifo minimize 
the resonant amplitudes in order to reduce roughness, noisenes§ 
and susceptibility to fatigue failure. The principal 
sources of dissipation of energy from the blade are : 

(a) Friction at the root and the hub , 

(b) Internal damping due to inelasticity of the material, and 

(c) Damping due to the presence of surrounding gas ( Aerodyna- 
mic Damping ) . 

The dissipation of energy at the root due to the 
relative motion between the blade and the rotor in mechanical 
attachments. At lower rotor speeds this energy dissipation 
may be appreciable, however, at very high rotor speeds 
the centrifugal force essentially tightens the blade 
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sufficiently to eliminate the dissipation. Some analysis of 
this type of damping has been reported in Refs. [190-192], 
Goodman and Slump [l9l] have pointed out that when blade 
root and hub are joined in press fit, the energy is dissipated 
by the microscopic slip at the interfaces. An experimental 
determination of the damping at the locking joints, when the 
turbine blades are vibrating, has been carried out by Eozolov 
[ 192 ]. By investigating the different types of locking 
joints-, the important parameters affecting damping have been 
determined. 

Although these investigations offer some understanding 
of the phenomenon of root damping, it is very much dependent 
on the type of the root, manufacturing tolerances, speed of 
rotation etc. Hence it is difficult to account for it in 
the analytical calculations. 

The material damping must often be relied upon, where 
structural slip and other types of external damping are not 
effective [193]. The material damping is caused by inelasticity 
of the material. The stress strain curve for loading and 
unloading are not generally identical but form a hyst crisis 
loop. The area bounded by the loop is a measure of the 
ehergy dissipated in each cycle. A brief resume of the 
damping capacity of the materials is given by lazan [194], 

Some investigation, of the effect of internal damping in 
rotating beams has been carried out by Morduchow [l95]* 
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An extensive study of the flexural and torsional 
vibrations of beams, with a special reference to turbine 
blades, taking into account the internal damping has 
been carried out by Pisarenko [196]. Baker et al. [l97] 
have considered the vibrations of cantilever beams in 
surrounding air* The damping arising from the air drag as 
well as from "the internal material hysteresis have been 
taken into account. 

Effect of internal damping on the transverse vibration 
of non-uniform cantilever beam is considered by Kerlin [198]. 
Grinin and Maybee [l99] has analyzed the effect of internal 
material damping on the stability of rotating systems. 

The aerodynamic damping is caused by the fact that the 
velocity acquired by a vibrating blade changes the incidence 

of the air stream. If the cascade characteristics are such 
that an increase of incidence causes an increase in the lift 
force on the blade, damping takes place. However, if an 
increase of incidence causes a decrease in the lift, it 
results in flutter, Pearson [200] has pointed out that in 
the unstalled region, the damping present is almost wholly 
aerodynamic, even when the excitation itself is almost 
wholly aerodynamic, e.g. from wakes, etc. He has derived 
an expression for the energy removed by the airstream from 
a vibrating blade per unit time. Cavaille et al.[20l] 
have presented a study of aerodynamic damping of turbomachine 
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blade vibration. 

Jones et al. [202] have described vibratihg beam 
dampers for reducing vibrations in gas turbine blades. 

Pai’in and Jones [203] present a preliminary investigation 
of a tuned damper utilizing a polymeric damping material 
encapsulation to withstand centrifugal load and erosion 
problems , 

2.19 MATERIAL 

The gas turbine imposes on materials the most severe 
condition known to date. Beside extreme mechanical strength 
at temperatures upto approx. 1000°G, corrosion resistance 
to a highly reactive environment must also be maintained. 

This controversial combination of properties has posed a 
Challange that has been successfully met by material 
scientist, has been done since Vf'orld War II, part iciilarly on 
Jet aircraft turbine. 

Improvement in the physical properties of materials 
have always been enthusiastically received by the aeronautical 
engineer. Basic knowledge in aerodynamics and structures 
has always been ahead of the metallurgical developments, and 
the designer can always make use of improved materials as soon 
as they become available. This situation is particular true 
at the present time. The lack of materials having the required 
physical properties at elevated temperatures seems to be the 
most important handicap in the development of supersonic 
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aircraft gas turbine, and progress in this direction could 
be seriously hampered if research in the field of high 
temperature materials .are not very actively purstied. 

The field of high temperature materials is by no means 
restricted to those materials having very high melting points. 

High temperature problems in materials exist over entire range 
of temperature . 

The report [204] gives the current state * of - the 
art of material usage in aircraft gas turbine engines# 
including auxiallary power units, amd discusses the trends 
in future .’materials such as composite, powder metullargy, 
controlled solidification, eutectic alloys, refractory 
metals and ceramics. Hew processing techniques, such as 
thermo -mechanical processing (TMP) and Gatorizing TM are 
also discussed. 

The tendencies manifesting themselves in the evalution 
of alloys which are currently used in the construction 
of turbines are reviewed [205]* The techniques of powder 
metullargy and oriented solidification should make it 
possible to make better use of Hi alloys. The perspectives 
of application of composite materials are examined. If the 
resistance to impact of these materials appears insufficient 
in case of moving compressor blades, their other mechanical 


characteristics permit a major lightening of other parts less 
exposed to shocks. Oriented eutectics, alloys of particular 
composition which take up a composite material structure 
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by unidirectional solidification present a set of properties 
which leads to expectation that some will be used for 
making turbine blades or vanes. Huff and Betz [206] also 
considered directionally solidified eutectic alloys and 
their application to turbine blades and gas turbine 
components. They conclude that lot of further investigation 
is necessary especially with respect to casting technology. 

Jahnke and Brnch [207] also contend that composite 
structures consisting of high strength fibres or plates 
in ductile materices with outstanding high temperature 
properties are achievable in directionally solidified 
eutectics. This new class of materials represent a major 
innovation in gas turbine blade technology. The advantage 
and limitations of two more'' promising eutectic system and 
the relationship of these properties to turbine 
blade design are discussed, 

Endres [208] has taken design examples of some 
typical high temperature parts such as vanes, blades 
and segments shielding hot surfaces and showed what 
conditions are encountered during service. 

The results of impact test [209] on large, solid 
fiber composite (graphite and boron/epoxy) fan blades 
for aircraft turbofan engine are discussed. It also showed 
that composite blades may be able to tolerate iceball 
and small bird impact, Habotova and Shipov[210] carried out 
a series of tests with geometrically similar steel and 
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duralumin rotor blades on a compressor stage model. It is 
shown that, with the same level of blade freq^uency differences, 
the rotor with steel blades has approximately 2.8 times 
hi^er critical pressure at which flutter arises than 
the rotor with duralumin blades. 

A Finite element method of rotating ceramic gas 
turbine blade consisting of an airfoil and root is discussed 
[ 211 ] . Three dimensional results using both isotropic and 
orthotropic material properties and two dimensional study 
of root considering the effect of friction, contact area 
and root geometry is included. 

Schallar, Visser and Lien [212] has presented a 
three-dimensional analysis of a ceramic rotating blade. 

The analysis reveals that under a centrifugal force load the 
maximum stress occurs in the root of the blade. The root, 
a critical area for this blade, is then analyzed in detail 
in an effort to reduce this stress. This paper is the 
first step in a feasibility study to ascertain if ceramics 
are applicable to rotating blades for heavy duty gas turbines. 

G-a 3 rbor [213] has presented wide range of creep 
resistant nickel - base alloys now available and developed 
primarily for use at elevated temperatures in gas turbine 
engine. 



CHAPTER III 




PORMULATIOH OP THE PROBLEM 

3.1 Introduction 

In this chapter the transfer matrix, approach is used 
to determine "Sb mode shapes and natural freq_uencies of a 
turbine blade. The low temperature side of the turbine can be 
idealized as a large aspect ratio cantilever beam tapered 
in both depth and width. Some of the predominant factors 
like rotation, shear deflection, rotary inertia, root flex - 
ibility, stagger angle etc., which affect the response 
of the blades are considered. ’State matrix' is derived 
for a small element of the beam and this is used in the 
determination of the elements of transfer matrix. The 
Rtinge-Kutta method is used for finding the transfer matrix, 
f’or the beam subjected to pure bending. 

3.2 Formulation of State Matrix for Piu'e Bending Mode 

(3.W 

The slope of the ce.ntre line of the beam is 

dx 

affected by the bending moment M and the shear force V. 

The angle between the perpendicular to the face and the 
centre line of the beam is caused by the shear force acting 
on the beam. Inspection of Fig, (3.2.1b) shows this angle 
to be eq.ual to relation between this angle 

and the shear force causing it is 

V(x) = SAg(x) ( H +t) 


(3.2.1) 
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where G-Ag(x)= &A(x)/K:g is the shear stiffness, being a 
shape factor depending on the shape of the cross-section, a(x) 
is the area of cross-section. 

The usual bending relation for a beam is 

M(x) = EI(x) (3.2.2) 

Now the rotating beam is subjected to distributed centrifugal 
force field. Considering an element dx dy dz with coordinates 
(x,y,z) subjected to components of the centrifugal force per 
unit volume along x,y and z axis are given by Ref. [174]. 
si (x+r) 

Ry = ^ ri (y .Cos^0 - z SinO Cos Q ) 

R- = PS'i ( -y Sin0 CosG + z Sin^0 ) 
where f is the mass per unit volume of the beam. 

Now these above components are integrated between limits 
y= -D/2 to D/2 and z=w-^ to w+^ to obtain the forces on the beam 
element of length dx situated at x from root. These are 
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dR =-p.(x) w SinGCosg dx 
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dR = p(x) w Sin% dx (3.2.3) 
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where iJ.(x) is mass per unit length. 


It is assumed that taper in both depth and width is linear, 
D(x) = Dq- ; Depth variation 

W(x) = » Width variation 

where D^, ^o’^hlS depths and widths of the beam at the 

fixed and free ends respectively. 

Force balance gives (Ref .fig. 5. 2. 2) - 
-R+(R+dR)+dF^ = 0 
dR = - dF 

X 

and P 

V-( V+dV) ~dF^ .Cos p(x)(4 w(x) .dx = 0 

For small values of v|/ , 

= -p(x)c:o^ w(x)- p(x)A^ w(x) Sin^^ 
or 

= -p(x) i<J +n^ Sin%).w(x) 

Moment balance gives 

M-(M+dM) - (R + dR).dw - dF^.|^ + (7+dV) .dx + 

X. c. 

[dF^+p(x)o»^ w(x).dx].|^ - 

p(x)ji^(x)CO^’^ (x) .dx = 0 
Reglecting second-order terms and arranging, 

dw 2^„^ ^..2 


(3.2.4) 


dM(x) 

dx 


V(x) - p(x)iy(x)6i -^(x) 


(3.2.5) 


where 


R(x) 


jb 

/ 


\^i^) rC (x+r) dx , 


the tension in the beam element due to rotation. 
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Rewriting these equations as 
dw(x) 


dx 


-1. i|,,(x) + l/GA (x). 7(z) 


= 1/Bl(x) . M(x) 
dx 


dM(x) 

dx 


= (l-R(x)/(JA„(x) ) . V(x) + [ R(x) - 


= -H(x) +xf S±n% ). w(x) 


'^(x) 

^ (x) 


n(x)iy(x)«or ] 


O? 


6 


( 3 . 2 . 6 ) 


where p. is mass per unit length, is frequency 

of vibration of rotating beam,f\is speed of rotation, 

r is radius of the disk, L is length of the beam and 
2 

i is radius of g 5 rration. 

t/ 

1 

When the state equations (3.2.6) are written in the 
matrix form ; 


-w(x) 

t^(x) 

M(x) 

7(x) 


d 

dx 


-w(x) 


0 

1 

0 

Q(x) 

'f(x) 

qBPH'-- 

0 

0 

P(x) 

0 : 

M(x) 


0 

RR(x) 0 

RM(x) 

V(x) 


AW(x) 

0 

0 

0 


(3.2.7) 


RM(x) = [ 1 + R(x) . Q(x) ] 


where 



57 


AM(x) 


I^N(x) 

= -^i(x) i^(x) cZ-fROO 

P(x) 

= l/EI(x) 

Q(x) 

=-l/GAg(x) 


or 

[zCx)] = [ B(x) ] [z(x)] (5^2*8) 

where [ z(x) ] is the state vector and [B(x)] is the 
matrix of variable coefficients of system of first-Kjrder 
differential equations (5.2,6). 

The above state vector and the state matrix would 
be used in the formulation of transfer matrix method which 
follows in the next section, 

3.3 THMSFER MATRIX 
A. Description 

The matrix which connects the properties ( the 
elements of the state vector) at 'n’ th section to that 
of the properties at 'n + 1' th is called a transfer 
matrix. The name 'transfer matrix' itself implies that 
to those of the properties of one section are related 
to another section through that matrix. This manner of 
description is eminently suited to describe the performance 
of a series of structures connected in tandem ; that is, 
output of the first is the an* put to second -jr, and so on. 

The transfer characteristics of such a Sicrios is found by 
multiplication of the transfer matrices of the individual 
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structures, i.e., local transfer matrices. 

Transfer matrix applies only to the case of one - 
dimensional structures, i.e., structures depending on 
just one spatial co-ordinate, and of certain two-dimensional 
structures where separation of spatial variables is 
possible. Consider a beam undergoing a harmonic motion, 
the state of motion at a given station on the beam may be 
described by a vector. 

e 

If excitation is about from the interval between stations 
a and b on the beam, then 


Z 


w 

i u 


= T 


b>— a 


(3.5.1) 


^b*— a transfer matrix associated with the structural 

element connecting station a to station b, a 4 X 4 matrix 


in this case . For the transfer of h dimensional state 
vectors ( k = even, for structural engineering), the 
transfer matrix is k X k. 

When excitations are introduced at a station, say 
at p, the value of the state vector becomes discontinuous 
at p. If p lies between a and b, equation (3.3.1) must 
be modified to accoiint for the discontinuity. Thus, 


^o ^p 


(3.3.2) 
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where is the input at p, o and IT are end stations. 

Since two components of and two of Zj^ are zeros 
(upon imposing boundary conditions, this is valid even when 
an end is, e.g., elastically constrained hy a spring, in 
which case the end station is considered to he immediately 
beyond the elastic spring, i.e, the transfer matrix, T must 
include the effect of elastic spring), equation (3.3.2) 
can be used to determine uniquely Z^ and Zj^. Knowing these 
we can return to equation (3.3.1) and compute;- the response 
at any station on the beam. 

The impedance, admittance, reoeptance and transfer 
matrices contain identical information, but in different 
arrangements and their interrelationships are given in 
Ref. [15]. 

B. Some Important Properties of Transfer Matrix 
(a) Cross ~ Symmetry 

If the structural element represented by a 

transfer matrix is syinmetrical about centre of element, 

by a suitable choice of the order and sign convention for 

the components of a state vector, it is always possible to 

obtain a transfer matrix symmetrical about its cross-diagonal 

(that is, the diagonal connecting "tij-,! * ^k-1,2 

ti - 1 . where t, . are the entries of the transfer matrix), 

1 »^ 

This cross -symmetic property can be expressed as 
^i,0 “ ^k+1-3, k-+l-i 


(A) 



(b) Inver s ion 
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From E(i.(3.3.1), 

= ^Wa (E) 

Tliis means "th.a'fe "the inverse of a "transfer mafrix 
is also a "transfer mafrix, and "bh.a"t "blis pair may be viewod 
as two propagation mechanisms in opposite directions 
provided by the same structural element. Furthermore, 
the existence of an inverse guarantees that a transfer 
matrix is non-singular. Dropping the subscripts of transfer 
matrices for brevity and denoting the elements of by 
^ijf then the inversion of a cross-symmetric transfer 
matrix merely amounts to change of signs for certain elements, 
namely 

^iO “ ( -1 ) ^ (C) 


For non-S 3 rmmetrical structural elements for which transfer 
matrix is no longer cross-symmetric, (C) becomes 



^k+l-j,k+l-i 


(D) 


( c) Value of Determinat# 


It follows from Eq.(C) or more generally Eq.(D) that 
I T I = j T”^ I = 1 (E) 

(d) Characteristic Equation and Eigenvalues 

By use of Eqns. (D) and (S), it can be shown that 
the characteristic equation of a transfer matrix I and its 
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inverse T ^ are identical and they can be written as 

+ + aKjd^- +1=0 

^ + 2K^3’‘-^ + 2K^3^ - 2Kj_3 + 1 =0 

Here k is even for structural engineering applications. 
Therefore T and T”^ have the same set of eigenvalues. Since 
for every eigenvalue of T there is a corresponding 

eigenvalue = l/o for T"^, we further conclude that 
the complete set of eigenvalues of T( or of T"^) must be 
pairwise reciprocal, 

(e) Helta - matrix of a Transfer Matrix 


TheA matrix of a transfer matrix also possesses 
property (a) stated earlier, i.e., the A matrix of a cross- 
S 3 raimetric transfer matrix is cross -symmetric. However, its 
inversion does not enjoy the same relationship as that in 
Eq. (D). Instead it satisfies. 


( T^)’^ = ( T”^ 

,hu,-l 


(G) 


so that I T I = I (T { = 1, that the characteristic 
equation of T"^ has a symmetrical form of the type of Eq,(F) 
and that the eigenvalues of T'^ also form reciprocal pairs. 


C. Frequency Equation and Mode Shapes 


The transfer matrix, by definition, is independent 
of boundary conditions and is a function of frequency , <»> , 
Upon substituting appropriate boundary conditions for the 
system into the state vectors of equation (3,3»1)» a 
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frequency determinant can be obtained. 

The systems now discussed are 'open' systems in 
the sense that there have been two sets of B.C's. If the 
system is closed, however, the boundary conditions are 
not so obvious. Suppose that we start at an arbitrary 
point of the closed system, where the state vector is 
z^. With the help of transfer matrices it is possible 
to proceed around the system until we return to the original 
point from which we started. In this way a relation of 
the form 


z = U U T. 
0 n n-1 


-U. 


TJ U 
2^1 


= T, z 
1 o 


- I ) -o = 0 

is obtained, so that the frequency condition is 

I Tl - ll= 0 

The frequency determinant is, therefore, of order n 
of the system and not of order n/2 , which is the case 
with 'open' systems. 

Having determined the natural frequencies, the 

coresponding normal mode shapes can also be obtained. This 

a 

is accomplished by relating one of the non zero components 
of the state vector at free end to 1.0 to provide a check 
on the mode shape computations. Then the state vectors 
are determined at other points through the system, by applying 
the transfer matrix. 
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D, Derivation of Transfer Matrix by the Runge-Kutta Method 
The Runge-Kutta method is well suited for finding 
the transfer matrix, but by and large it is necessary to have 
the services of an electronic computer. To compute the 
transfer matrix relating the state vectors and we 

define the interval between o and RS into a number of small 
intervals. Let us suppose that between sections n and n+1 
the interval size is h (Rig. 3.3.1). 

The Eq. (3.2.8) is rewritten as 

[z(x)] = [B(x)] [z(x)] (3.3.3) 

The Range - Kutta method gives the resiilt 


Section® ^ 2 


Local 

Transfer 

Matrix 


U. 


z 


U. 


n n +1 
i t 


RS-1 RS 


h 1 


!^nt 


I * 
I f 


Ri 


R3 


o -^1 2 


^n ^n+1 


®RS-1 %S 


Eig.(3.3.1) The subdivision of the interval between o and RS, 


=n+l = S I ^ 


where 


(3.3.4) 

K. 


% = hB(x^+|) ) 


K 

Kp = hB(x +|)(z^+^) K 3 = hB(x^+h) (z^+K^) 

(3.3.5) 

B(x ), B(x +h), and 3(x +^) are, respectively the 

values of the matrix B(x) at sections n and n+1 and at 
the section midway between n and n+ 1 . 

Substituting expressions (3*3.5) into E(i. (3.3*4). 
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Thus the relation between z and z ; 

n+1 n 

^n+l = I ^ I [ + 2 ) + B(x^+h) ] + ^ 

B(x^ + |).B(Xj^) + B(x^+h) .B(xj^+|) + B^(x^+|)] 

3 

+j^[ •B(Xj^) 4- B(x^+h).B^(x^+|-) ] 

4- |^B(x^+h).B2(x^+|).B(x^) 

(3.3.6) 

where the contents of the braces represent, the local 
transfer matrix linking the state vector and • 

By repeating the process for all the intervals between o and 
FS, the overall transfer matrix T relating z^ and Zjjg can 
be found by multiplication of the braces of Eq., (3.3»6) in 
the proper sequence. In Big. (3.3.1) it is clearly shown 
that the overall transfer matrix can be obtained by Successive 
multiplication of these local transfer matrices. 

From Big. (3.3.1), 

^1 ^1^0 » ^2 ^2^1 ' •♦^FS ^FS^FS-1 

(3.3.7) 



By backward 
Therefore, 


where , 


substitution Zj^g is related to z^ . 

I^IfS ~ I^IfsI^IfS-1*’** "* I^U 


z 


T 


FS 


= T z 


(3.3.8) 


= . U 


FS 


^Ifs-1 


uigluli 
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5 *4 Determination of Natural Frequencies a.Tid .ghapPH 



Writing Eq. ( 3 . 

.3.1) in 

expanded 

form 

9 

-w 


'^11 ^12 

%3 

"14 



t 


'^21 ^22 

^23 

"24 


t 

M 


^31 ^32 

"33 

"34 


M 

V 


^^41 ^42 

"43 

^44, 


7 


NS 


o 

(3.4-1) 


where the index 'NS' corresponds to the free end and 'o' 
corresponds to the fixed end of the team. 

Imposing the boundary conditions that is shear and 
bending moment would vanish at free and (NS), and deflection 
and bending slope would vanish at the fixed end ('o'), 

' 0 

0 

M 

Y 

(3.4.2) 

The above equation leads to the following condition : 


1 

-w 

■df 


T 

0 

= 

0 

NS - 


"11 

'^12 

*13 

*14 

"21 

^22 

"23 

*24 

’^31 

'^32 

"33 

*34 

1*41 

^42 

"43 

*44j 




(3.4.3) 


Non— trivial solution of the above equations exist only 
when the corresponding determinant vanishes . 



Therefore, 


T 

T 


T 

35 34 

T 

43 44 


0 


(3.4.4) 


The values of ’oa * for which Eq..(3.4.4), called frequency 
equation, is satisfied are the natural frequencies of 
the system. These natural frequencies are determined 
by a linear interpolation procedure. 


Mode Shanes : 

Eq. (3*4.2) gives. 


- Wjjg - T^^ Vq 



To, 

23 0 


+ T 


24 


V 

o 


(3.4.5) 


Assuming deflection at free end as unity, using 
Eqs. (3.4.3) and (3.4.5) the following relations can be 
obtained . 


State Vector at Free End : 


■■ ^NS 

1.0 

II 

CO 

dajT; 

%s = 

0 


0 

State Vector at Fixed End 


" ^ 24^33 ^'^^^ 13^34 ” ^ 14 ^ 33 ^ 


(3.4.6) 
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X = 0 

~'^ 33 ^^^ 13^34 


T14T33) 


T T} ) 
U 33^ 


(3.4.7) 


This initial state vector, preaiultiplied by successive 


local transfer matrices, would lead to state vectors at the 
ends of successive segments. So, the values of deflection 
at each segment end would lead to the mode shape. This 
procedure is applied to all the natural freq.uencies to get 
the crossponding mode shapes. The requirement that - w^jg = 1.0 
at the tip provides a check an the mode shape computations. 


3.5 Root Flexibility Effect : 

In a turbomachinery the blades do not form an 
integral part of the disk, but are rather loosely mounted 
on the rotor disk. The vibrations of the blade, therefore, 
continue to some extent into the mounting. This has an 
effect similar to that of increased blade length and 
results in a lower frequency than the one calculated on the 
fixed e«.nd assumption. The root flexibility varies with the 
speed of rotation, since the centrifugal forces tend to tight- 
en the blade in the disk and increase the fixidity. 

The root flexibility, generally makes the root a 
non-deflecting end, but not a non-rotating end ; therefore, 
imposing boimdary conditions 

“ns = ''sS - ° ’'0 = ° • 'I'o = ' 


in Eq. (3.4.I) gives 
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"l 2 *14 


0 



*21 *22 *23 *24 


M/K 

0 


*31 *32 *33 *34 


M 


NS 

.*41 *42 *43 *44. 


V 


(3.5.1) 

where K is rolling stiffness factor. Obviously 
very high value of K leads to fixed end condition, whereas 
very small value lead to simply supported end. 

The above equation leads to the following condition ; 
(T32A + T33) = 0 

(T^g/E + 1^3) = 0 


Non-trivial solution of the above equations exist only 
when the corresponding determinaiit vanishes. 

Therefore, 


( T 
( T 


32 ^^ + ) 

^34 

42 /K *43 ’ 

'^44 


0 ( 3 . 5 . 3 ) 


The values of ' for which the Eq. (3*5.3), called frequency 
equation, is satisfied are the natural frequencies of the 
system* 

Mode Shanes : 

Bq. ( 3 . 5 JL) gives, 

- Wjjg = ( 1^2/^ ■*' ^13 ^ ^0 ^14 
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Assuming deflection at free end as unity, using Eqs. (3.5.2) 
and (3. 5 *4) the following relations can he obtained. 

State Vector at Free End ; 




1.0 


[ ( ) + K ( 

0 

0 (3.5.5) 


State Vector at Fixed End : 

- w 

o 

to 

«o 

7 

o 

where 

T = = x^2 •^34 ” -^14 -^32 ^ ^ -^13 ^34 ” -^14 ■*■33 

Effects of root flexibility can easily be taken into 
account by taking various values of K, such as 10 ^ 10 , 

10 etc. covering the whole range from fixed root 
condition to simply supported root condition. 


0 


*34/-' 

E *34^ 

- ( Ij2 + K Tjj ) / T 


(3.5.6) 


f m m 


\ I XT' ^ rn 


m \ 
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OHiPTBR 17 


4.1 Formulation of State Matrix for Combined 
Bending and Torsion 


The deflections at the centre of mass and the shear 
centre are connected by the ecjuation. 

v(z) = w(x) + e . P (x) ( 4 . 1 . 1 ) 

where v(x) = deflection at the centre of mass of beam 
cross-section at distance 'x' from the root, ■w(x) = 
deflection at shear center, p(x) = torsional deflection, 
and e = distance between center of mass and shear center. 
Considering force equilibrium for stagger angle 

*0 ’ , we have 

~V(x) dx + ti(x)<iO^ v(x) dx + MSI? v(x)Sin^a dx = 0 

( 4 . 1 * 2 ) 

Substituting Bq. (4.1.1) into Bq. (4.1.2), 

J^(x) =z (CO^ +O^Sin^i9 ) •w(x)-n(x) ( ckD^+ri^Sih^0 ) e.p(x) 

(4.1.5) 

Equilibrium of torsional moment gives 

“t~ ( |?^T M.e +qi(x)4oV(x).p(x) = 0 

Jr 

(4.1.4) 


or 

(x) = p.(x) et^^. w(x) + |i(x)CO^ [ i^(x) + e^] . ^=t(x) 
dx y 

(4.1.4) 


Consider the following equation : 

^l(x) [ iy(x) + e^ .] = ^(x) Ip(x)/A(x) 

(4.1.5) 
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where = po3.a:.r moment of inertia and i.(x) = area of 

cross-section. 

Substituting Eg., (4-. 1.5) in Sq. (4.1.4), we get 

|_iljj(x) _ e. w(x) + p.(x)ck^ Ip(x) / a(x) .p(x) 

( 4 . 1 . 6 ) 

Flexural moment equilibrium leads to the following equation, 

dM(x) = l/( x) - p.(x) i^(x)Cc:>^ " "4^(3:) - R.(x) . 
dx J I 

[ w(x) + e .p (x) ] 

(4.1.7) 

Slope '^1^' is ginresn by 

= -q.(x) + l/SAg(x) . 7(x) (4.1.8) 

Simple Bending iCenaent theory gives the following relation, 

= l/Bl(x) . M(x) (4.1.9) 

dx 

Torsional deflection 'p' is related to the torsional 
moment 'M,p' by the following equation, 

= l/dJ(x).M^(x) (4.1.10) 

where ' J' is mass moment of inertia about centre of mass. 
Rewriting coupled bending-torsion equations, 

||(^> = -v|.(x) + l/SA^Cx) . V(3C) 

— = A/BI(x) . j9(x) 

dx 

= [ R(x) - p(x)iy(x)€4>^ ] ..sj.(x) 4-V(z). 

[ 1-R(x) /OA„(3c) 3-R(x)e •^,p(3c) 

*" "S 
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= -ii(x) (Co ~ r(x) icJ- +Sl%^)e,^U) 

I (x) 

^ p.(x)cc? e. w(x) + ti(x)6a^ ’AU) * 


=l/&J(x).M,^ 


(4.1.11) 


Writing the sta,te equations (4.1.11) in the matrix form: 


1 

^-w(x> 


0 

1 

0 

Q(x) 

0 

0 


-w(x) 


'^(x) 


0 

0 

P(x) 

0 

0 

0 


(x) 


M(x) 


0 

RE(x) 

0 

RM(x) 

AT(x) 

0 


M(x) 

dx 


•f 










7(x) 


AM(x) 

0 

0 

0 

0 

AU(x) 


V(x) 


M,3,(x) 


AS(x) 

0 

0 

00 

0 

AV(x) 


,M,3,(x) 


P(x) 


0 

0 

0 

00 

AW(x) 

0 


P(x) 


( 4 . 1 . 12 ) 

where L 

R(x) =j \x{x) ri (x+r) dx 

mix) = |i(x) ieo^ +ri^ Si]fe) 

iHir(x) = -!J-(x) iy(x)6»^ + R(x) 

RM(x) = [ 1 + R(x) .Q (x) ] 

P(x) = 1 /EI(x) ' 

Q(x) =-l/(xAg(x) 

AS(x) = - \xix) ^ e 
AT(x) = - R(x) . 5 j(£)' 

AU(x) = -AM(x)- . e 

a7(x) = ia(x) £0^ /a(x) 

= 5J(x) 



or 
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l^j-EzCx)] = [c(x)] . lizCx)] (4.1.13) 

where [z(x)] is the state vector and [G(x)] is the state 
matrix of variable coefficients of first-order coupled 
bending-torsion differential equations (4.1.11). 

The Runge-Kutta Method is again used to obtain the 
transfer matrix, which will be of the order 6X6, 

4.2 Determination of Natural frequencies and Mode Shanes 


Writing the transfer matrix in the expanded form, 





C\J 

H 

^13 

*14 

^15 

*16 





*21 

^22 

^23 

*24 

^25 

*26 



M 


'^31 

T52 

'^33 

*54 

^35 

*56 


M 

V 

— ■ 



'^43 

*44 

^45 

*46 


V 

Hji 


Si 

®52 

^53 

*54 

^55 

*56 





.’'61 

^62 

^63 

*64 

"65 

*66. 




NS 


( 4 . 2 . 1 ) 

Taking the root flexibility into account and imposing 
the boundary conditions 

“hs = T^IIS = %s = ° V^o 


in Eq. (4.2.1) 
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.J 


fp 

■^11 

Ti2 

^13 

®14 

*15 

*16‘ 


'O' 



^21 

T22 

^23 

*24 

*25 

*26 


M/K 

0 


*^31 

^32 

^33 

*34 

*35 

*36 


M 

0 


^1 

^42 

^43 

*44 

*45 

*46 


Y 

0 


^51 

"52 

T^3 

*54 

*55 

*56 


% 

P 


?61 

^62 

^63 

*64 

*65 

*66. 


0 


( 4 . 2 . 2 ) 


0 

NS 

The above equation leads to the followijog condition : 

( 1^2/^ + T35 > + T34 + T35 = 0 

( T42/K + T43 ) + T45 = 0 

( + T53 ) + T54 7 ^ + T55 = 0 (4.2.3) 

Non-trivial solution of the above equations exist only when 
the • corresponding determinant vanishes. 

Therefore , 

( T32/K + ’^33) 

( T^g/K + 1^3 ) T 44 

( T^g/K + T53 ) ’^54 ''55 

The values of ■ t» ' for which the M. ( 4 . 2 . 4 ). called 
frequency equation, is satisfied are the natural frequencies 
of the coupled tending- torsion system. 

Mode Shanes : 


"34 


35 

'45 


= 0 ( 4 . 2 . 4 ) 


Bq. (4.2.2) gives. 
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- w 


US = ^ *13 ’ “0 *14 '^0 ^ ®15 ®T 

'•’hS = ^ *22^^ *25 ' “° *24 '^o *25 “t 

^HS = ^ * 62 /*^ * *63 ’ "0 + *64 * *65 “l 


0 

(4.2.5) 

Assvunins deflection at free and as unity, using Bqs. (4.2.3) 
and (4.2.5) the following relations can he obtained. 


fpo.tor at Fre^ 

En( 3 . : 

“ ^NS 

= 1.0 



= A2/B5 

. T55 



• ^54 

^FS 

= 0 


%S 

= 0 


■^FS 

= 0 


^FS 

= Ag/B^ . 

. T55 


B^/Bg . 

■ ^54 


•25 


+ * ^55*^64 


■65 


( A5/35 + 

(4.2.6) 


State Vef-tnr at 'Fixed. Enj • 


-Wo -"0 


= I/B5.K >55 




1/B, 


55 


= V®6 ^55 




= - (A5/B5 + 
= 0 


“^ 54 ^ 



where 


•^L 

=: 

( 1 j^ 2 /K Ii 3 ) 

^2 

= 

( *22/^ + 5^23 ) 

h 

sc 

^ ^ *^33 ^ 

H 

- 

( ^42/"^ ^ "^43 > 


= 

( T52/K + T53 ) 


= 

( ^62/^ ^ "^63 ^ 


=: 

-^1 ^55 “ S ^15 



A5 T35 - A3 T55 

S3 


'^14 ^55 “ ^15 '^54 

S4 

— 

'^34 ^55 ' ^35 "^54 

S5 

= 

+ B2 . ^3/^4 

S6 

ss 

% B4 + 
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5.1 RESULTS AND COHCLUSIOIS 

In this study transfer matrix approach is used 
to obtain the natural frequencies and mode shapes of 
large aspect ratio blades. The effect of the parameuers 
like taper, root flexibility are also studied. The effect 
of shear deformation and rotatory inertia on the natural 
frequencies, which are of secondary importance are also 
invest igated . 

5 . 2 . Snecif icat ion s of the Blade 


Length of blade 

L 

= 50.0 cm 

Depth of blade at root 

d 

= 5.0 cm 

Width of blade at root 

b 

= 1.3 cm 

Radius of the disk at the 
blade root 

r 

=12.5 cm 

Stagger angle 

0 

= 1.3963 

Young* s modulus 

E 

0 2 

= 21 X 10° gm/ cm 

Modulus of rigidity 

G- 

0 2 

= 8.3 X 10 gm/em 

Mass per unit length 

P 

= 0-0516 gm/cm 

Moment of inertia 

A 

— about z axis 


= 13.54 cm^ 

• — about y axis 


= 0.9151 cm^ 
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5 2.1 T^ptermination of the Position of Shear Centre 

The position of shear centre e, from the centroid 
of the beam of singly symmetric rectangular tubular crosssection 
as shown in figure below, is given by^ 
d(B2 - Bj) 2b(B^ /t2- yti) 


e = 


2(B2 + B^) 

^2 




a 


Where B^, B^ — lumped areas at the corner of beam X-section 


t ''2 


— wall thicknesses 



The shear centre positions for t o.4 cm, 
b = 1.3 cm and varying values of t 2 are listed belo 


S.No. 


'2 cm 


cm 


— Off -Set distance 
in o/o age of Chord 



0.10 
o;a3 
0.50 
0.55 
0.60 
0 .65 
0.70 
0.75 
0.80 


0.0000 
0.0773 
0.1444 
0.2037 
0.2504 
0.3055 
0 . 3499 
0.3911 
0.4295 


0 o/o 

2 o/o 

3 o/o 

4 o/o 

5 o/o 

6 o/o 

7 o/o 

8 o/o 

9 o/o 


Taking + valuea for computation J- 

no of-Foct -. of off-set distances* e' on ® 


can evaluate 

frequeneies 
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5 . 3 Humer ical Co mpiit at io n 

Computer programs are written for the evaluation of 
natural frequencies for pure bending and coupled bending- 
torsion mode based on the formulation given in Ch. Ill and 
IV. Fortran IV language is used for the programs which 
are run on IBM 7044 computer. Flow chart is given in 
Appendix A., 

Effect of various values of double taper ratio 
on the variation of frequencies of first four ' bending 
mode^, taking different number of segments are obtained for 
three values of rotational speeds (5000,10000 and 15000 RPM) . 
These are presented in Table 1 a,b and c. Bending 
frequencies of uniform beam for different number of segments 
are also obtained and presented in Table 2, for various 

Values of rotational’, speeds (5000,10000 and 15000 RPM). 

\ 

Convergence studios have been made and plotted as in Figs. 
5.1.1 a,b both for non-uniform and uniform beam in bending 
mode at 5000 RPM. It is observed that the rate of 
convergence is faster and decreases with increasing modes. 

Bending, modes (first four) of uniform beam at 3000 
RPM are presented in Table 3 and plotted in Fig. 5.1. 2. 

First bending and bending-torsion mode (with 3 o/o 
off-set distance ' e' ) of uniform beam at 3000 RPM are 
given in Table 4. The effect of o'ff-sst distance is shown 


in Fig. 5.1.3 
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It is observed from the flgvxe that the displacements 
are close to the pure bending mode for distances close to 
the support. However, as we move away from the support, 
they start diverging with the maximum value occuring at 

the tip. 

Banding mode (seoond>f uniform Beam for different 


rotational speeds (3000,6000, and 10000 EPM) are given in 
fable 5. These results are plotted in Pig. 5.1.4. ™hioh 
shows that as the rotational speed inoreases say from 
3000 BPM to 10000 EPM, the beam region between 20 o/o to 
SO o/o is affeoted much more than that near the tip region. 

The reason being that tip deflection is constrained to 

unity for computational purposes. 

Coupled bending-torsion mode of vibration of uniform 

beam is obtained for five sets of off-set distance (e). 
fhese are presented in Table 7, for 3000 EPM and plotted xn 

Fig. 5.1*5. 

. • /n-p f-iY'qt four bending natural 

The variation of first lour _ ^ 

•th the chano-e in the stagger angle ( ) of e 

frequencies with the cna » 

t'qM p 6 An increase of 0 from 
root of the blade is shown m Table . ^ 

0° to 90° decreases the fundamental freiuenoy. by abou 

r>-ipq orenot significantly affected. 
14.5 o/o. The other frequencies orenot s . 

These are plotted in Fig. 5 •l*5a. 

4- -PI AYlbilities on the frequencies 

The effect of root flexibilities 

of first four bending modes of 
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ratio = 0.5) as well as that of uniform beam are presented 
in Tables 8 and 9 by taking various values of rolling 
stiffness factor K (eg. 10^®, 10^^, 10^ and 10“®), for various 
values of rotational speeds (0,5000 and 10000 jElPM) . These 
results are plotted in Pig. 5*1.6 a and b, which show that as 
rolling stiffness factor is increased, there is a marked 
increase in the natural frequencies ■w[ith the increase in 
the rotational speed. The results obtained for various values 
of K in the limit coincides with those obtained for a 
cantilever beam. 

Bending frequencies of non-uniform (double taper ratio 
= 0.5) as well as uniform beam for various values of rotational 
speeds (from 0 BPM to 10000 RPM) are also obtained to show 
their effect and presented in Table 10. 

The variation of frequency parameter ratio with depth 
and width taper are obtained for various values of rotational 
speeds (0,5000,10000 and 15000 SPM) and presented in Table 11. 

The above results for the first three modes are plotted 
in Pigs. 5*1.7 a,b and c and are compared with the results 
of Rao and Carnegie [9l]* It can be seen that they are in 
good agreement , account ing for shear deformation and rotary 
inertia effects, for non-rotating case. 

For the non-rotating case, for a given width taper 
the frequency parameter ratio for the fundamental mode 

decreases as -depth taper increases. For the second and third 

■ * 



81 


modes, the frequency parameter ratio shows a marked increase 
with the depth taper. Also for a given depth taper the 
frequency parameter ratio decreases as the width taper 
increases for all the three modes. 

The effect of rotational speed is also shown 

in Fig- 5.1.7. 
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RBCOiMEJTDATIONS 

The following are some of the points on which 
further investigations ace recommended to he carried 
out . 

1) Use of transfer matrix approach for obtaining 
natural frequencies of twisted, tapered and rotating 
blades with and without root flexibility. 

2) The effect of disk elasticity. ' 

3) The effect of camber on the natural frequencie 

4) Extension to small aspect ratio blades. 

5) Study of Flutter characteristics. 
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1“ 1 o nti “blip V 811*1.3,1^ ion of f npi^iisncy wiuii 

Segments fSS various values of rotational speeds. 


no . 


iioMJTMiTpnHM BR« (BBSBIIia MOBB ) 


/I = 5000 RPM 


double taper ratio 


0.25 



— — HO 

6 

"Of SEOMERTS 

10 12 

20 

25 

30 

953.410 

2160.173 

4046.542 

6710.714 

952.50 

2156.650 

4039.729 

6700.366 

>48.141 

2138.207 

4017 .409 

6695.980 

-945.286 

2131.133 

4009.431 

6716.408 

, 930.795 

2098.588 

4028.257 

5613.935 


DOUBLE TAPER RATIO 


0.50 



710.252 

2119.824 

4424.219 

7980.466 

716.722 

2131.146 

4487.227 

7951.831 

717.950 

2154.814 

4484.998 

7916.789 

719.793 

2140.497 

4488.785 

7883.333 

720.174 720.178 

2141.769 2141.822 
4490.707 4493.411- 
7882.302 7882.352 

DOUBLE taper RATIO 

= 

0.75 



575.030 

2175.341 

5012,378 

8934.273 

577.522 

2176.039 

A947.745 

8964.821 

578.002 

2176.953 

4940.388 

9016.160 

578.726 

2178.792 

4935.905 

8965 .729 

578.876 

2179.241 

4936.079 

8960.384 



table 1 (a) 
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NOI-UniPOiM BEM ( BESDINff MODE ) 
-Q- = 10000 RPM 


DOUBLE TAPER RATIO =0.25 


Uo. OF SEOMEUTS 

6 10 12 20 25 


1712.350 1711.765 1705.645 1702.528 1749.352 

3525.684 3522.574 3491.852 3480.493 3508.060 

5925.634 5903.824 5851.852 5463.291 5233.306 

9525.684 9569.825 9601.852 9611.528 9643.306 


DOUBLE TAPER RATIO = 0.50 


1269.823 

3358.786 

6243.401 

8962.151 

1271.016 

3357.414 

6179.518 

9822.376 

1272-827 

3362.191 

6177.232 

9897.232 

1275.778 

3370.789 

6175.055 

9898.009 

1276.404 

3573.193 

6181.053 

9866.386 

DOUBLE TAPER RATIO 

=r 

0.75 


998.840 

3337.857 

6604.857 
1 0526.667 

999 .707 
3335.167 
6530.366 
10931.259 

1000.437 

3336.145 

6500.394 

10881.932 

1001.594 

3338.697 

6516.417 

10846.348 

1001.848 

3339.243 

6515.465 

10835 .710 


TABLE 1 (b) 



85 


HOH-UNIPOK'd BSm (BEIDIITG MODE ) 




n 

15,000 RPM 


DOUBLE 

TAPER RATIO 

0.25 


Mo 

. OF SEUMEMTS 



6 

10 

12 

20 

25 

2512.143 
4325.645 
6725.978 
1 0325.978 

2511.233 

4322.042 

6703,744 

10369.744 

2505.189 

4291.644 

6651.144 

10401.989 

2503.008 

4280.727 

6263.958 

10411.445 

2567.907 

4308.060 

6033.306 

10443.308 

DOUBLE taper RATIO 


0.50 


1854.442 
4655.742 
8421.740 
1 6645.740 

1857.407 

4677.407 

8539.407 
15460,409 

1958.000 

4693.750 

8600,417 

14960.417 

1861.429 

4632.929 

8632.928 

13823.928 

1872.667 

4797.071 

8890.571 

11142.904 

DOUBLE TAPER RATIO 

= ; 

0.75 


1413.684 
4652.449 
8479.480 
1 8624.480 

1434.187 

4638.187 
8445.245 

18556.245 

1436.966 

4651.739 

8439.239 

17202.863 

1439.583 

4661.196 

8427.862 

13359.239 

1437.734 

4659.347 

8492.661 

12955.181 


TABLE 1 (c) 
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Variation of frequency with no. of segments for various 

values of rotational speeds. 

UNIFORM BUM (BENDINU MODE) 

a = 5000 RPM 
No. Of Segments = 25 


6 

10 

12 

20 

25 

■ 30 

478.238 
2238.627 
5431.433 
1 0045 .142 

478.024 

2232.101 

5393.274 

10116.017 

478.008 

2231.498 

5382.948 

10053.230 

477.996 

2230.984 

5372.551 

9986.684 

477.995 

2230.934 

95 371. 541 
9980.766 

477.996 

2230.935 

5371.540 

9980.765 

0 = 10000 RPM 

798.656 
3345.969 
6,969.754 
1 1692.021 

797.518 
3340.247 
6387.774 
11874 .751 

797.433 

3339.418 

6876.328 

11809.415 

797.386 

3338.714 

6865.050 

11750.466 

797.388 

3338;7I9: 

6864.^30 

.11744.502 


O. = 15000 RPM 

1132.893 

4616.241 

8810.241 
13049.241 

1124.049 

4609.635 

8779.655 

13083.635 

1124.979 

4607 .820 
3768.156 
14205.656 

1123.354 

4606.975 

8759.975 
19509.976 

1125 .152 

4606.569 

8759.815 

21284.816 



TABLE 2 
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Bending Mode Shapes of Uniform Beam at 5000 RPM 

MODE SHAPES 


Int erval I 

End 

II 

III 

1 

0.0000 

0.0000 

0.0000 

2 

0.0037 

-0.0156 

0.0413 

3 

0.0143 

-0.0568 

0.1434 

4 

0,0307 

-0.1163 

0.2764 

5 

0.5233 

-0.2639 

0,4156 

6 

0.0785 

-0.3402 

0.5347 

7 

0.1086 

-0.4109 

0.6185 

8 

0.1422 

-0.4715 

0.6527 

9 

0.1788 

-0.5179 

0.6305 

10 

0.2182 

-0.5469 

0.5518 

11 

0.2578 

-0.5561 

0.4226 

12 

0.3034 

-0.5438 

0.2548 

13 

0.3486 

-0.5094 

0.0644 

14 

0.3953 

-0.4529 

- 0.1296 

15 

0 .4432 

-0.3854 

- 0.3076 

16 

0.4921 

-0.3785 

-0.4508 

17 

0.5417 

-0.2418 

-0.5438 

18 

0.5920 

-0.0865 

- 0.5753 

19 

0.6426 

++0 .0834 

M ) . 5400 

20 

0.6936 

0.2616 

-0 .4383 

21 

0 .7447 

0.4445 

-O .2762 

22 

0.7959 

0,6297 

-0.0648 

23 

0.8470 

0.8151 

+0.1819 

24. 

0.8991 

0.8542 

+ 0 .4494 

25 

0.9491 

0.9005 

+ 0.7247 

26 

1.0000 

1.0000 

1.0000 


IV 


0.0000 
-0.0819 
-0.2634 
-0.4676 
-0.S291 
-0.7004 
-0.6571 
-0.5005 
- 0.2561 
+ 0 .0326 
+ 0.3122 
0.5305 
0.6464 
0.6379 
0.5073 
0.2801 
0.0008 
-0.2753 
-0.4927 
- 0.6058 
-0.5876 
-6.4366 
-0.1617 
+0.1935 
+0.5918 
1.0000 


table r3- 
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Effect of Off-Set Distance ( 5o/o ) on the first 
Bending Mode Shapes of Uniform Beam at 3000 RPM 

MODS SHAPES 


lETBRVAli PURE BENDING BENDING-TORSION 

END 


1 

0.0000 

2 

0.0049 

3 

0.0193 

4 

0.0428 

5 

0.0749 

6 

0.1156 

7 

0.1629 

8 

0.2172 

9 

0.2774 

10 

0.3422 

11 

0.4106 

12 

0.4814 

13 

0.5530 

14 

0.6342 

15 

0.6936 

16 

0.7593 

17 

0.8202 

18 

0.8748 

19 

0.9218 

20 

0.9601 

21 

0.9889 

22 

1.0080 

23 

1.0173 

24 

1.0177 

25 

1.0110 

26 

1.0000 


0.0000 

0.0051 

0.0226 

0.0476 

0.0871 

0.1402 

0.2082 

0.2929 

0.39&0 

0.5200 

0.6676 

0.8427 

1.0494 

1.2940 

1.5838 

1.9285 

2.3441 

2.8382 

3.4419 

4.1817 

5.0677 

6.2394 

7.6807 

9.5171 

11.8804 

1.0000 


TABLE 4 
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Yariation of Bending Mode Shapes (Second Mode) with 
Rotational Speeds. 


MODE SHAPES 


Interval 3000 RPM 6000 RPM 10000 

End 


1 

0.0000 

0.0000 

2 

-0-0156 

-0.0178 

3 

-0.0568 

-0.0630 

4 

-0 .1163 

-0.1269 

5 

-0.2639 

-0.2024 

6 

-0.3402 

-0.2834 

7 

-0.4109 

-0.3647 

8 

-0.4715 

-0 .4419 

9 

-0.5179 

-0.5107 

10 

-0.5469 

-0.5678 

11 

-0.5561 

-0.6099 

12 

-0.5438 

-0 .6346 

13 

-0.5094 

-0.6395 

14 

-0.4592 

-0.6232 

15 

-0.3854 

-0.5847 

16 

-0.3785 

-0.5233 

17 

-0.2418 

-0.4396 

18 

-0.0865 

-0.3342 

19 

+0.0800 

-0.2089 

20 

0.2616 

-0.0659 

21 

0.4445 

0.0922 

22 

0.6297 

0.2662 

23 

0.8151 

0.4410 

24 

0.8542 

0.6253 

25 

0.9005 

0.8122 

26 

1.0000 

1.0000 


0.0000 

-0.0199 

-0.0671 

- 0.1301 

-0.2013 

-0.2753 

-0.3478 

-0.4155 

-0.4753 

-0.5248 

-0.5615 

-0.5833 

- 0.5882 

-0.5746 

-0.5412 

-0.4869 

-0.4113. 

-0.3145 

-0.1972 

- 0.0608 

0.0923 

0.2595 

0.4374 

0.6222 

0.8106 

1.0000 


TABLE 5 
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Effect of Sta,i;ger Angle variation at 3000 RPM 

NON - UNrPORM 3ENDIN0 MODE) 

DOUBLE TAPER RATIO =0.5 
NO. of Segments ••= 25 


o 

o 


6 ^ 60° 

a=9o° 

605.577 

584.557 

540.786 

517.353 

1789.348 

1781.412 

1768.475 

17.61.445. 

4039.192 

4036.148 

4030.027 

4026.949 

7376.347 

7375.046 

7371.758 

7369.035 


UNIFORM BEAM 

(BENDING MODE) 


477.067 

1932.958 

4995.562 

9560.126 

450.500 

1926.594 

4993.061 

9559.249 

392.018 

1913.728 

4988.118 

9556.608 

359.120 

1907.243 

4985.626 

9554.290 


TABLE 


6 
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Effect of off-set difjtance on the variation of frequency 
COUPLED BENDING - TORSION MODE 


UNIFORM BEAM 
/I = 3000 RDM 


No. of Se, aments 


= 25 


0 o/o 


OFF-SET distance IN PERCENTAGE OF CHORD 
3 o/o 5 o/o 7 o/o 9 o/° 


453.568 

2389.931 

6323.158 

12226.977 


453.583 

2390.000 

6307.653 

12175.150 


546.718 

2390.326 

6279.447 

12118.014 


547.800 

2390.369 

6310.715 

12243.515 


547.825 
2400.526 
6331.715 
12189 .549 


table 7 
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Effect of ioot Flexibility on frequency for various values 
of rotational speeds. 


NON-UIIPORM BEAM (BEBBIEG MODE) 


DOUBLE TAPER RATIO 

a 


= 0,5 

= 10000 RPM 


ROLLING STIFFNESS FACTOR , K 

j_q3B ^q20 ^q4 

10~® 

1276.404 

1276.404 

1173.944 

1173.969 

3575.193 

3373.193 

3141.766 

3142.256 

6181.053 

6181.053 

5726.058 

5725.801 

9866.386 

9866.386 

9184.712 

9139.434 




n = 

5000 RPM 

720.174 

2141.769 

4490.707 

7882.302 

720.174 

2141.769 

4490.707 

7882.302 

594.285 

1817.652 

3872.718 

6968.774 

594.279 

1817.645 

3872.686 

6967.693 



a = 

0 RPM 

356.129 

1505.129 

3738.624 

7065.125 

356.129 

1505.129 

5738.624 

7063.125 

236.514 
1016.669 
2977 .451 

6035.462 

236 .528 

1016.646 

2977.432 

6 035. 404 


TABLE 8 
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Effect of rioot Flexibility on frequency for various 
values of rotational speeds. 


OTIFORM BEiM (BBIDIEG MODE ) 



a 

= 10000 RPM 



ROLL INI STIFFNESS FACTOR 

, K 

1058. 

10^° 

10^ 

10"® 

797.388 

3338.719 

6864.230 

11744.502 

797.388 

3338,719 

6864.230 

11744.502 

663.999 

3084'.058 

6289.727 

1Q785.215 

664.003 

3084.087 

6289.638 

10784.680 


Cl 

= 5000 RPM 


477.995 

2230.934 

5371.542 

9980.766 

477.995 

2230.934 

5371.542 

9980.766 

332.387 

1869.379 

4607.372 

8821.133 

332.381 

1869 .369 

4607.360 

8821.199 


/I 

= 0 RPM 


271.177 

1698.810 

4753.240 

9306.508 

271.177 

1698.810 

4753.240 

9306.508 

185.395 

1189.734 

3854.822 

8037.185 

185.401 

1189.734 

3854.822 

8037 .185 


table 9 
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Variation of frequency with the rotational speeds 

NON-OTIPOBM BEAM (BENDIFG- MODE) 
DOUBLE flPER RATIO =0.50 
No. OP SEOMEMT = 25 


0 RPM 


2000 


4000 


6000 8000^ 


10000 


356.129 

1505.129 
3733.624 
7063 •IS^ 


500 

1700 

4100 

7400 


750 

2250 

4350 

7650 


1050 1250 
2550 3150 
4950 5550 
8250 9150 


1276.404 

3373.193 

6181.053 
9866 .886 



UNIPORM BEAM 

(bending mode ) 


271.177 

1698.810 

4753.240 

9306.508 

316.264 
1894 . 87 4 
4858.239 

9618. 032 

413.398 

2055.796 

5158.660 

9743.013 

540 ^ 700- 

2427.218 3100 
6619. 409 6400 

10280.648 11200 

797.388 

3338.719 

6864.230 

11744.502 


table 10 
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l/’a.ria’t ion of fr©(5[u©nc'/ pfi.i'aiUG'tsi’ nal/io wifli dspili and 
width taper. 


fl 


Co 


o 


= 0 RPM 

= 271.177 CPS 


Mode 

Width 

Depth 

Taper ,6 



Paper, 

Y 0.25 


0.50 

0.75 


0.25 

2.733 


2.402 

2.240 

1 

0.50 

2.001 


1.726 

1.587 


0.75 

1.654 


1.383 

1.264 


0.25 

27.415 


35.131 

42.927 

2 

0.50 

23.993 


30.796 

37.657 


0.75 

22.244 


28.566 

34.919 


0.25 

136 .703 


201.466 

268.464 

3 

0.50 

123.133 


189.547 

253.221 


0.75 

123.754 


183.578 

245.734 



TABLE 11 

a 





= 5000 RPM 




COq = 477.995 CPS 



0.25 

3.79 


2.61 

, 1.765 

1 

0.50 

3.52 


2.28 

1.515 


0.75 

3.35 


2.18 

1.435 


0.25 

19.2 


22.3 

22.7 

2 

0.50 

17.2 


20.998 

21.5 


0.75 

16.1 


20.291 

20.8 

3 

0.25 

0.50 

0.75 

70.55 

66.50 

63.95 


91.50 

88.50 
87.00 

108.195 

105.762 

103.200 


P ABIE 11 h 
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variation of frequency parameter ratio vith depth end 
width taper. 

a = 10000 RPM 

= 797.388 CPS 

^0 


Mode 

Width 

1'o*nPT • Y 

Depth Taper 

n-28 

,5 

.0^50 

0.75.^ 

« • 

1 

0.25 

0.50 

0.75 

4.797 

4.797 

4.797 

2.89 

2.55 

2.39 

1.87 

1.665 

1.57 

2 

0.25 

0.50 

0,75 

18.3 

18.3 

18.3 

18.45 

17.82 

17.72 

18.72 

17.90 

17.50 

3 

0.25 

0.50 

0.75 

43 

43 

43 

61.70 

60 

59.25 

68.81 

67.32 

66 .50 


sable 11 c 


A 

CC> 0 


15000 EPM 
1125.152 CPS 


1 

0.25 

0.50 

0.75 

5,200 

5 .21 

5.25 

3.01 

2.77 

2.74 

2 

0,25 

0.50 

0.75 

14.70 

14.30 

14.10 

18.35 

18.20 

18.10 


- 

— — * 

62.00 

61.50 

61.40 

3 

0.25 

0.50 

0.75 

28.70 
27.60 
. 27.10 


1.71 

1.65 

1.63 


17 .65 
17.32 

17.20 


58.00 
57 .38 
57.10 


table 11 ^ 
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FIG. 2.13.1— COMPRESSOR MAP SHOWING BOUNDARIES FOR 
FOUR TYPES OF FLUTTER. 
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FIG, 3.2.2 -FORCE BALANCE DIAGRAM. 
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FIG. 5.1-1- b- VARIATION OF FREQUENCY WITH NUMBER OF SEGMEh 
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